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ABSTRACT 


J) An  automated  cloud  analysis  program  vas  developed  and 

/" 5~cf»  h  A.  [JV-f..  Pr*<- 

established  on  the  SP ADS 'computer  system  at  "the  Naval  Envi- 
ronaental  Prediction  Research  Facility  (NEPRF) .  The  pro* 

v  fTjCpsf a-^iongJ'V  Opt  rein*  *-L  £\vir^r^P.n+s.l  fh  ft )  y 

graa  evaluates  GOES^  visual  and  infrared  satellite  iaagery 
simultaneously.  The  analysis  aethod  produces  inforaation  on 
cloud  zypes,  cloud  amount,  precipitation  intensity ,  and 
cloud  top  height  and  teaperature  through  use  of  threshold 
tests  of  radiance,  texture,  and  teaperature.  A  review  of 
current  work  on  the  evaluation  of  satellite  inforaation  by 
computer  and  by  manual  analysis  is  included. 

A  maritime  region  460  X  463  nautical  miles  in  size  was 
selected  for  test  analysis.  The  satellite  imagery  was  manu¬ 
ally  evaluated  and  compared  to  the  computer  generated  out¬ 
put.  Reasonably  good  patterns  of  cloud  types,  precipitation 
and  cloud  amount  were  produced  by  the  computer,  although 
further  testing  and  verification  is  needed. 
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i.  imsuusxxsi 


The  personnel  operating  Naval  ships,  aircraft  and  weapon 
systems  need  near  real-time  sab  synoptic  and  sesoscale  mari- 
tiae  weather  inforaation  to  operate  effectively  and  safely. 
Current  conventional  observational  data  (froa  surface  ,  raw- 
insonde,  and  aircraft  observations)  do  not  adequately  cover 
the  aaritiae  region  cr  provide  for  aesoscale  analysis  of 
weather  phenoaena.  Specifically  needed  are  detection  and 
aonitoring  of  any  adverse  weather  conditions.  This  need 
pertains  to  phenoaena  ranging  froa  aesoscale  (10-100  ka)  to 
synoptic  scale  (100-1000  ka)  which  are  poorly  resolved  by 
conventional  observations.  Satellite  observations  can  solve 
this  nroblea  through  the  acquisition  of  global  high  density 
data  which  are  available  regardless  of  local  influences  of 
geography,  surface  conditions,  or  local  politics. 

Another  problea  is  the  need  to  provide  an  analyzed  prod¬ 
uct  to  the  user  in  a  short  (1-3  hour)  time  frame.  Manual 
analysis  of  satellite  data  is  time  consuming  and  tedious, 
and  it  does  not  provide  the  information  in  a  real-time 
fashion.  An  interactive  computer  system,  on  the  other  hand, 
can  collect,  analyze,  and  produce  the  required  products 
within  the  required  time  limit. 
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This  research  effort  concentrated  on  the  development  of 
techniques  to  specifycloud  types,  cloud  amount,  cloud 
height,  and  precipitation  intensity  with  an  emphasis  on  the 
aaritiae  regions  froa  satellite  images.  Information  on 
critical  weather  parameters  such  as  low  visibilities,  ceil¬ 
ings,  precipitation  presence,  and  intensity  can  be  derived 
froa  these  cloud  characteristics.  This  work  uses  Geosta¬ 
tionary  Operational  Environmental  Satellite  (GOES)  Visual- 
Infrared  Spin  Scanned  Sadioaeter  (VISSB)  digital  satellite 
data  from  visual  and  infrared  channels  which  have  a  one-half 
hour  temporal  and  0. 5  to  4  nautical  miles  (n  mi.)  spatial 
resolution.  The  data  were  received  and  processed  at  the 
Naval  Environmental  Prediction  Research  Facility  (NEPRF)  in 
Monterey,  California  using  the  Satellite  Data  Processing  and 
Display  System  (SPADS).  This  interactive  computer  system 
was  developed  at  NEPRF  for  use  by  the  operational  Navy  to 
give  real-time  information  about  weather  phenomena.  One  sys¬ 
tem  is  now  operational  at  the  Naval  Eastern  Oceanography 
Center  (NEOC)  at  Norfolk,  Virginia. 

The  objective  of  this  thesis  is  to  establish  and  evalu¬ 
ate  a  computer  program  on  the  SPADS  to  analyze  GOES  images 
to  derive  cloud  types,  precipitation  intensity,  cloud 
height,  and  cloud  amounts.  To  accomplish  this,  an  in-depth 


study  of  previous  research  and  analysis  techniques  was 
initiated,  and,  froa  this  study,  useful  algorithms  and 
techniques  were  gleaned.  These  algorithms  and  techniques 
were  combined  to  produce  the  desired  computer  program.  Pre¬ 
liminary  tests  of  the  computer  program  were  completed  on  the 
SPADS  system  at  HEPS? . 

Chapter  II  begins  with  a  review  of  previous  studies. 
This  includes  discussions  of  four  techniques  of  cloud  type 
classification;  (1)  cloud  index  (Harris  and  Barrett,  1975, 
1978),  (2)  spectral  analysis  (Liljas,  1981a,  1981b),  (3) 

3-dimensional  nephanalysis  (Fya,  1978)  ,  and  (4)  two-dimen¬ 
sional  histogram  (Platt,  1 981)  .  A  discussion  of  the  techni¬ 
ques  for  determining  cloud  height  and  cloud  amounts  used  by 
Reynolds  and  Yonder  Haar  (1977)  and  others,  and  of  proce¬ 
dures  for  identifying  precipitation  areas  and  their  associ¬ 
ated  intensities  is  also  covered.  Chapter  III  describes  the 
set  of  algorithms  selected  for  implementation  on  SPADS. 
Chapter  IV  presents  an  outline  of  the  test  procedures  used 
and  reports  the  results  of  these  tests  (including  detected 
problem  areas).  Chapter  7  summarizes  the  research  efforts 
and  problems  and  provides  suggestions  for  future  changes  and 
additional  research. 
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A.  INTRODUCTION 

Satellite  visual  (Vis)  and  infrared  (IR)  imagery  can 
provide  a  vide  variety  of  weather  information,  including 
estimates  of  cloud  type,  cloud  height,  cloud  aaount,  and 
precipitation  intensity.  In  fact,  a  satellite  nephanalysis 
can  produce  a  spatial  suaaary  of  these  iaportant  cloud  fea¬ 
tures.  An  analyst  exaaines  the  visual  brightness  to  derive 

estiaat9s  of  thickness  of  the  cloud  and  the  cloud  texture  to 

• 

differentiate  between  cuaulus  and  stratus  type  clouds.  By 
using  the  IR  iaagery ,  cloud  top  heights  can  be  estimated; 
the  brighter  the  cloud  iaage,  the  colder;  therefore,  the 
higher  it  is.  With  coabined  visual  and  infrared,  one  has  a 
greater  ability  to  pick  out  the  cloud  types  and 
furtheraore,  to  estiaate  the  probability  of  precipitation. 

One  of  the  aain  problems  encountered  in  the  1960 *s  and 
early  1970  's  was  the  lack  of  standardization  of  rules  gov¬ 
erning  the  nephanalysis  construction.  Each  country,  and 
soaetiaes  each  analyst,  has  used  subjective  methods  for  pro¬ 
ducing  cloud  analyses.  Some  of  the  specific  prcbleas  aris¬ 
ing  froa  lack  of  standardization,  as  discussed  in  Harris  and 


Barrett  (1975)  ,  have  bees:  no  minimum  size  is  established  as 
reference-  for  inclusion  or  exclusion  from  analysis;  cirras 
clouds  are  rarely  identified  in  aiddle  latitudes; 
significant  cloud  areas  are  subjectively  delineated  at  the 
discretion  of  the  analyst;  only  four  cloud  categories  are 
assigned  (which  consisted  of  unequal  percents  of  areas)  ;  and 
no  standard  nephana  lysis  construction  rules  have  been 
published. 

With  these  problems  in  Bind,  Harris  and  Barrett  (1975, 
1978),  Barrett  and  Martin  (198  1),  Platt  (1981),  Fye  and 
Logan  (1977)  and  Lil  jas  (1981a)  ,  to  name  a  few,  have  been 
working  toward  the  goal  of  producing  a  recommended  standard 
for  objective  nephana lysis  of  satellite  VIS  and  IS  imagery 
through  manual  and  computer  analysis.  The  following  sections 
review  these  representative  research  efforts  in  this 
endeavor.  The  discussions  are  divided  into  three  areas  of 
concern  to  this  thesis;  cloud  typing,  cloud  height  and 
amounts,  and  precipitation. 

3.  CLOUD  TYPING 

Basically  four  distinct  approaches  to  effective  cloud 
typing  have  been  proposed  in  the  literature.  They  are:  (1) 
cloud  index  (Harris  and  Barrett,  1975,  1978;  Barrett  and 


Harris, 1977)  ,  (2)  the  (J.S.  Air  ?orce  3DHBPH  program  (Fye  and 
Logan,  1977;  Fye,  1978),  (3)  spectral  cross-correlation 

(Liljas,  1981a,  1981b),  and  (4)  two-dimensional  histogras 

(Platt,  1 98 1)  .  These  eethods  are  reviewed  below, 

i.  Qaa£  lAisi 

Harris  and  Barrett  (197  5)  recognized  the  need  for 
standardization  and  fa:  increased  information  content.  As  a 
result,  they  developed  a  cloud  indexing  systen  which  pro¬ 
vided  specific  guidance  and  rules  for  analyzing  satellite 
visual  ieagery  manually.  For  clarity  and  understanding, 
they  proposed  the  analysis  be  composed  of  three  separate 
layers  of  analyses;  two  descriptive  and  one  interpretive. 
Fig.  1  is  the  flow  chart  of  the  stages  in  the  construction 
of  this  nephanalysis.  The  authors  developed  a  key  (Fig.  2) 
which  includes  details  on  the  following  cloud  features;  (1) 
percent  of  cloud,  (2)  cloud  type  (  e.  g.  cumulus),  (3)  cloud 
structure  (including  size,  shape,  and  pattern),  and  (4) 
height  of  cloud  tops.  The  key  also  gives  interpretive 
weather  phenomena. 

The  first  analysis  layer  contains  information  on 
geography  and  cloud  cover.  The  cloud  cover  is  divided  into 
five  equal  percent  categories  with  a  minimum  size  reported 


?h«  second  descriptive  layer  includes  cloud  type 
and  cloud  structure.  They  are  classified  is  the  following 
order:  eajor  cloud  systees#  definite  boundaries  and  finally 
indefinite  boundaries.  The  analyst  uses  cloud  brightness 

and  texture  to  discern  six  types  of  cloud3  (see  Fig.  2). 

i 

The  final  analysis  layer  produced  is  the  interpreta¬ 
tion  of  the  cloud  features  using  the  information  from  the 
first  two  layers.  The  types  of  cloud  feaxures  considered 
iaportant  are  displayed  in  Fig.  2.  All  three  analysis  layers 
coabined  were  designed  to  be  coapleted  within  one  hour  by  a 
skilled  analyst. 

Barrett  and  Harris  (1977)  extended  this  procedure  to 
include  infrared  imagery  for  the  purpose  of  supplementing 
the  visual  analysis.  This  added  capabilties  for  nighttime 
and  high  latitude  analyses.  They  mention  three  possibilities 
for  analyzing  IR  imagery;  (1)  use  longstanding  internation¬ 
ally  accepted  codes  and  symbols,  (2)  follow  the  visual  pro¬ 
cedures  previously  outlined,  or  (3)  use  a  new  procedure 
which  recognizes  and  represents  the  special  properties  of 
the  IS  images.  The  first  two  methods  would  be  directly  com¬ 
parable  to  the  visual  analysis.  The  third  approach  would 
offer  additional  information  to  allow  identification  of 
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Figure  2.  VIS  Key  to  be  Osed  in  Manual  Naphanalysis  (from 
Harris  and  Barrett*  1975). 
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cloud  height  through  cloud  brightness  evaluation  and  a  more 
specific  identification  of  cloud  texture. 

The  IB  analysis  technique  identifies  the  following 
cloud  features  (see  Fig.  3):  (1)  percent  of  cloud  cover 

(same  as  Vis  technique),  (2)  brightness  of  cloud,  (3)  cloud 
texture,  and  (4)  size  of  clouds  (same  as  VIS)  ,  plus  inter¬ 
pretive  categories  with  eighteen  of  them  the  same  as  VIS 
(heights  and  three  cloud  descriptions  are  not  included). 
Barrett  and  Harris  (19  77)  felt  that  any  of  the  three  methods 
could  produce  an  analysis  in  a  single  hour  if  done  by  a 


trained  analyst 


q ed 


Figure  3.  IH  Key  to  be  Use4  in  Kanual  Nephanalysis  (from 
Barrett  ana  Harris,  1  977). 

After  establishing  a  manual  method  for  both  7IS  and 
IH  satellite  nephanalysis,  Harris  and  Barrett  (1978)  devel¬ 
oped  an  autoaated  objective  nephanalysis  through  a  cloud 
recognition  scheme.  They  used  an  approach  which  examined 
quantitative  cloud  brightness  and  texture  within  a  small 
area  of  the  whole  picture.  The  flow  chart  in  Fig.  4  out¬ 
lines  their  approach.  Cloud  brightness  was  assessed  in  each 


subarray  by  suaaing  the  brightnasj  counts  of  only  cloud  pix¬ 
els  and  then  di elding  by  the  number  of  cloud  pixels  (2q.  1 

and  2,  Table  I).  Cloud  aaount  was  calculated  in  each  subar¬ 
ray  by  suaaing  the  nuaber  of  cloud  pixels  and  then  diriding 
by  the  total  nuaber  of  pixels  per  subarray.  This  nuaber  was 
then  multiplied  by  100  to  obtain  a  percant  of  cloud  corer 
(Sq .  3  in  Table  1)  . 

Their  assessaents  of  cloud  texture  U3ed  statistical, 
instead  of  structural.  Measures  of  local  rariation  in  image 
density.  One  measure  they  employed  was  the  standard  deria- 
tion  of  each  subarray  (5x5  pixel  array)  of  density  ralues, 
as  calculated  with  Eq.  4  (Table  I). 

Their  other  criterion  was  rector  dispersion,  which 
required  that  the  density  ralues  be  treated  as  a  set  of 
adjacent  triangular  planes.  The  dispersion  in  three  dimen¬ 
sions  of  the  normal  tc  these  planes  gave  another  measure  of 
texture.  Pig.  5  shews  a  graphical  representation  of  this 
method.  Eg.  5,  6,  and  7  in  Table  I  were  used  to  calculate 
these  estimates  of  texture. 

Once  all  these  ralues  were  calculated,  boundaries 
between  cloud  types  had  to  be  established.  Three  categories 
were  selected;  (1)  sheet  or  layered  clouds  (stratiform). 
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(2)  cellular /tower  clouds  (cumulus)  and  (3)  broken/aottled 
cloud  (stratocuaulus  or  aized)  (Harris  and  Barrett,  1978) .  k 
three-diaensional  decision  space  was  set  up  using  two  tex¬ 
ture  eeasure aents  and  one  brightness  value.  Those  areas  that 
had  a  brightness  above  the  threshold  were  included  in  fur¬ 
ther  analyses;  those  below  were  not.  Fig.  6  shows  the  deci¬ 
sion  space  based  on  the  two  texture  measurements.  The 
discriainant  lines  were  developed  through  analysis  of  a 
training  set  of  Defense  Military  Satellite  Program  (DHSP) 
saaple  brightness  and  texture  aeasureaents  of  each  cloud 
category . 

In  assessing  the  accuracy  of  this  method,  Harris  and 
Barrett  (1978)  found  several  problea  areas.  For  axaaple, 
towering  cumulus  in  post  frontal  areas  (especially)  were  not 
identified  because  of  the  field  of  view  (FOV)  size  and  snow 
areas  were  included  as  clouds.  But  the  overall  accuracy  was 
still  good  with  a  greater  than  72%  correct  classification. 
Harris  and  Barrett  (1978)  felt  that  several  procedures  would 
improve  the  output  of  the  model.  Among  these  were  increasing 
the  resolution  of  the  data  by  decreasing  the  size  of  scan¬ 
ning  spot,  decreasing  the  subarray  size,  normalizing  the 
image  brightness  for  changes  in  sun  viewing  angle,  and 
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increasing  the  texture  parameters  to  perait  identification 


Figure  5.  Graph  of  Vector  Dispersion  Technique.  (al 

represents  a  smooth  surface  (stratiform  clouds) 
ana  (b)  represents  a  rough  surface  (cumulus 
clouds).  Normal  vectors  for  each  case  are 
included  to  the  left  of  the  diagrams  (from  Harris 
ana  Barrett,  1978). 

data  from  aircraft,  surface,  and  upper  air  observations.  \ 
cloud  analysis  is  produced  in  a  real-time  fashion  regardless 
of  ice,  snow,  desert,  or  other  geographical  anomalies  (Fye 
and  Logan,  1977)  .  The  satellites  provide  high  density  data 
over  the  entire  world,  thus  permitting  a  global  analysis  of 
clouds.  The  horizontal  resolution  of  the  output  is  25  n  mi. 

The  data  array  has  a  pixel  resolution  of  3  n  mi. 
The  satellite  data  are  first  formatted,  rectified,  mapped, 
and  stored  on  a  computer  storage  device  before  running  the 
3DNEPH  program.  The  rectification  corrects  for  the  effects 
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Alogrithss.  for  the  OtjectiTe+He^^|&lysis  (froe  Harris  and 
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t  -  density  threshold  defining  the 
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A  -  average  cloud  amount 
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-  standard  deviation  of  each  sub¬ 
array  denalty  values 

-  aean  of  the  subarray  density  values 


n^,  a^)  -  direction  cosines  of  the  ith  plane 


(1,  n) 


(6)  R2  -  ait)2  +  (In^2  +  (Int)2 


-  nax  likelihood  of  true  polar 
vector  (X,  u,  y) 


(7)  k 
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8  ■  the  sanple  size 

k  •  estimate  of  the  texture  of  the  surface 

where: 


k-1  •  smooth  surface  (stratiform) 
k~0  •  uneven  surface  (cellular  cumulus) 
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of  liab  darkening  before  storage.  Liab  darkening  is  the 
effect  of  a  decrease  in  brightness  as  the  edges  or»  "limbs" 
of  the  earth  scan  are  approached.  The  data  are  not  normal- 
ized  with  respect  to  light  (VIS)  or  longwave  radiation  (IB) 
(Fye  and  Logan,  1977).  Each  visual  and  infrared  pixel  is 
represented  by  a  single  grayshale  value  in  the  range  froa  1 

to  63.  In  the  visual,  1  is  dark  and  63  is  white.  Infrared 

pixels  represent  temperature  in  degrees  Kelvin  (K)  from  210 
K  for  a  value  of  1  to  310  K  for  63;  each  grayshade  repre¬ 
sents  a  1.6  degree  change  (Fye,  1  978).  The  preliminary  pro¬ 

cessor  calculates  the  average  grayshade  value  per  25  t  ai. 
square  using  Eg.  1  (Table  II)  and  then  uses  this  to  get  the 
variability  within  the  square  using  Eq.  2  (Table  II) 
(Fye, 1978).  Both  these  calculations  are  performed  sepa¬ 
rately  for  IR  and  VIS  images.  The  VIS  and  IR  grayshade  val¬ 
ues  and  variabilities  are  used  to  identify  the  various  cloud 
types.  Fiq.  7  and  8  represent  the  cloud  typing  algorithms 
used.  They  were  constructed  empirically  by  comparing  many 
grayshades  and  variabilities  to  corresponding  satellite 
imagery  and  surface  data  (Fye, 1978).  The  operational  algor¬ 
ithms  were  derived  from  the  figure  and  then  modified  to 
obtain  optimum  results.  One  obvious  problem  with  electing 
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to  use  the  sodel  (giv®  the  available  docuee station)  is  the 
lack  of  graduated  values  on  the  axes  of  these  two  figures. 
Only  relative  values  of  vara/cold  or  light/dark  are  included 
in  the  report.  Also  lacking  in  the  Fye(1978)  report  are  any 
details  of  the  current  algorithas  used  foe  cloud  typing. 


TABLE  II 

3DNEPH  Statistical  Equations  for  the  Average  Grayshade  and 

Variance  (from  Pye,  1978) 


N 

(1)  G  •  1/H  G  ,  G>0,  N>0 

i-1 


G  -  average  gray  shade 

»  individual  3  nmi  area 
grid  gray  shade 

N  *  total  number  of  G^  per 
25  nmi  grid  space1 


N 

(2)  V-l/sJ]  |Gt  -  G|  ,  G>0,  N>0 
i-1 


V  -  variability  within  the 
25  nmi  grid  space 


3.  SBSslx&l 

The  third  classification  aethod  is  the  spectral 
cross-correlation  method  developed  by  Erik  Liljas  (1981a) 
for  use  on  an  interactive  computer.  Liljas  used  TIROS-N  and 
National  Oceanic  and  Atmospheric  Administraion  (NO AA)  6 
satellite  data  from  three  spectral  channels;  channel  1 
(VIS)  ,  channel  2  (near-IR)  and  channel  (IR)  .  During  the 


summer  time,  different  cloud  types,  land  and  water  surfaces 
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7.  The  Decision  Space  of  the  3DNEPH  Program.  IE 
and  VIS  gravshades  are  used  to  determine  the 
cloud  type  (from  Fye,  197  8). 
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were  analyzed  to  detect  their  unique  spectral  signatures. 
The  ra suits  were  used  to  build  a  three  diaensional 

parallelepipedon  classification.  Pig.  9. 

Two  satellite  types  were  used  far  the  test  of  the 
nodel;  polar  orbiting  (TIBOS-N  and  HOAA-6)  and  geostationary 
(Heteosat)  .  The  satellite  data  used  for  the  test  of  the 
aodel  were  geoeetrically  corrected,  but  they  did  not  need  to 
be  norealized  because  the  aid -day  pass  of  the  satellite  was 
used  (approximate  sun  elevation  of  45  degrees). 

Noraalization  corrects  for  the  radiance  variation  caused  by 
the  sun  location  (elevation).  Also  the  aaount  of  data  was 
reduced  froa  1,024  digital  levels  to  256  levels,  which  was 
considered  adequate  resolution  for  the  study.  Data  from 
channel  2  were  priaarily  used  to  separate  land  and  water 
areas.  Channels  1  and  4  were  the  principal  cloud  type 
identifiers. 

Coaparison  of  the  spectral  cross-correlation  method 
with  synoptic  observations  showed  that  the  computer  method 
gave  correct  indication  of  cloui  types,  and  also  an  improved 
cloud  division.  As  shown  in  Pig.  10  (see  Table  III  for  sym¬ 
bol  definitions  of  cloud  types) ,  several  cloud  types  are  in 


the  same  boxes; 


Liljas  recommends  the  use  of  texture 
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TABLE  III 


10  (from  Liljas, 


Main  and  Cloud  Typos: 

1.  Cumulonimbus  * 

a.  Scots  cloud  vtch  high  cop 

a£  Squall  cloud  with  scmctarad  showers 

2.  Mlabostratus  b 

b.  Large  vercical  thickness 
bj  Rather  low  topside 

2.  Cirrostratus  c 

c1  Dense  cirroscratus 
cf  Cirrus 

’cf  Thin  cirrus  over  water 
Dense  altoscracua 
di  Thin  altostracus  over  water 
e|  Thin  aleoscracus 

4.  Cumulus  conges tus  f 

a.  Dense  altocumulus  * 

f?  Large  piled  up  cumulus 

.  f j  Rather  mall  and  flat  ctaaulua 

5.  Stratocumulua  h 

h.  Dense  stracocustulus 
hi  Ordinary 

ff  Slightly  piled  up  cumulus 
J  with  dear  areas  in  between 

i.  Very  dense  haze/stratus 
ij  Dense  haze/stratus 

6.  Haze/Stratus  i 

i-  Ordinary  haze/stratus 
g  Cumulis  humilis 

7.  Land  1 

j  Haze  over  water 

1.  Planting  season  spring  or  autumn 

l£  Warm  green  season 

S.  Water  k 

It.  Cold 
k£  Warm 
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*•  iMordiagnsiQuai  qiataqtfti 

Another  spectral  method  vas  developed  by  Platt 
(1981)  using  a  two-dimensional  histograa  representation  of 
cloud  types.  He  selected  satellite  data  with  typical  cloud 
fields  to  develop  a  histograa  decision  space  representation 
of  cloud  types.  Platt  concentrated  on  evaluating  cloud  sys- 
teas  instead  of  isolated  clouds  and  therefore  his  results 
are  biased  toward  synoptic  features. 

The  measured  visual  radiances  were  compared  with  the 
direct  radiance  froa  the  solar  dish  on  the  satellite  to  pro¬ 
duce  a  calculated  isotropic  albedo  per  2.5  km.  These  albe- 
does  are  expressed  as  a  fraction  from  0  to  1.  The  infrared 
radiances  are  converted  to  blackbody  brightness  temperatures 
ranging  from  193.3  to  301.5  K.  A  histogram  was  then 
constructed  by  sorting  these  data  into  983  different  bins  of 
temperature  and  albedo.  Each  area  evaluated  was  approxi¬ 
mately  500  X  500  km,  with  roughly  40,000  values  per  histo¬ 


gram. 


The  total  bin  counts  were  plotted  on  a 


two-dimensional  grid  and  then  contoured  by  same  bin  count 


clouds.  Fig.  12  vas  compiled  to  give  an  idealized  depiction 
of  the  location  of  clcud  types  by  albedo  and  temperature. 

Platt  (1981)  discussed  several  factors  that  tend  to 
reduce  the  usefulness  of  Fig.  12.  One  cloud  problem  arises 
when  there  are  either  breaks  in  a  cloud  deck,  or  partially 
filled  field  of  view  (FOV)  in  the  satellite  view  area. 
Additionally,  clouds  of  the  sama  optical  depth,  but  exhibit¬ 
ing  breaks  will  have  a  different  albedo  measurement  than  an 
unbroken  cloud  of  the  same  optical  depth  (stratocumulus  in 
Fig.  14  demonstrates  this) .  Another  factor  that  changes  the 
albedo  and  temperature  values  is  variation  in  cloud  top 
height.  All  of  these  factors  era  ate  ambiguity  in  determining 
the  cloud  type,  amount,  optical  depth,  and  cloud  top 
temperature.  As  a  result,  the  author  recommends  that  addi¬ 
tional  work  should  be  done  with  these  results  before  using 
the  two- dimensio  ral  histogram  approach  to  identify  clouds. 
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Figure  11.  Two-dimensional  Histogram  Plot  of  Albedo  Versus 
Temperature  (representing  cumulus  clouds)  (from 
Platt,  1961)  . 
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Figure  12.  Ideal  Histograms  of  Cloud  Types  (from  Platt, 


C.  CLOUD  HEIGHT  &HD  AHOUHT 

Three  aethods  of  deriving  cloud  height  and  aaount  have 
been  proposed  in  current  research.  Tvo  aethods  are  an  exten¬ 
sion  of  those  reviewed  in  the  cloud  typing  section;  cloud 
index  (Harris  and  Barrett,  1975,  1978;  Barrett  and  Harris, 

1977)  and  the  U.S.  Air  Force  3DNEPH  (Pye  and  Logan,  1977; 
Fye,  1978).  In  addition,  a  new  bispectral  technique 
(Reynolds  and  Vonder  Haar,  1977)  is  discussed. 

1-  XH&2 

In  two  separate  studies,  3arrett  and  Harris  (1977; 
Harris  and  Barrett, 1 975)  developed  a  recommended  standard 
cloud  height  and  aaount  deterai nation  procedure  for  aanual 
analysis  of  satellite  visual  and  infrared  images.  The  vis¬ 
ual  method  combined  brightness,  cloud  type,  and  synoptic 
situation  to  determine  qualitative  estimates  of  cloud  height 
(low,  middle,  or  high  level).  Later,  they  developed  a  sys¬ 
tem  for  identifying  cloud  levels  through  the  degree  of  cloud 
brightness  in  th«  IB  image.  IR  cloud  information  specifi¬ 
cally  lends  itself  to  this  task  since  the  radiation 
temperature  of  the  cloud  is  a  function  of  the  cloud-top 
altitude  (Barrett  and  Harris,  19  77). 


Some  problems  arose  with  cirrus  and  cloud  eleaents 
smaller  than  the  field  of  view.  In  this  instance,  satellite 
aeasured  temperatures  can  appear  warmer  than  actual  cloud 
temperature,  because  the  satellite  is  also  receiving  the 
radiation  temperature  of  the  surface  or  lower  clouds.  This 
problem  can  be  minimized  by  using  18  and  VIS  together  to 
identify  cirrus  clouds  and  by  using  higher  resolution  data  , 
g.q.,  data  from  the  Advanced  Very  High  Resolution  Radiometer 
(AVHRR)  on  NOAA  satellites. 

Fig.  3  is  the  key  for  the  IR  nephanalysis  and 
includes  the  cloud  brightness  steps  and  Fig.  13  is  an  exam¬ 
ple  of  how  these  are  applied.  This  analysis  can  be  inter¬ 
preted  in  terms  of  qualitative  heights,  low  (dim)  to  high 
(very  bright)  (Barrett  and  ti arris,  1977)  .  cloud  amount 
determination  categories  are  also  delineated  in  Fig.  2  for 
VIS  and  Fig.  3  for  IR.  It  is  done  in  categories  of  size  of 
cloud  element,  e.g.  ,  100-150  km.  Harris  and  Barrett's 

(1978)  objective  automated  method  does  not  specifically 
address  cloud  height,  but  does  calculate  an  average  cloud 
amount  for  each  subarray  size  (5x5  pixels)  using  Eg.  3 


(Table  I) 
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The  U.S.  Air  Force  3DHEPH  program's  main  objective 
is  to  produce  information  on  cloud  height  and  cloud  aaount. 
It  uses  a  standard  vertical  grid  of  15  layers,  which  range 
fron  the  surface  to  100  aillibars  (ab|  .  The  IR  iaage  is 
used  to  deternine  height  and  ela aent  cloud  amounts,  and  the 
visual  and  infrared  are  used  together  to  determine  total 
cloud  aaount.  This  inforaation  is  than  combined  with  sur¬ 
face,  aircraft  and  rawinsonde  data  to  determine  the  final 
cloud  height  and  cloud  aaount  analyses. 

The  preliminary  processor,  as  described  in  the  cloud 
type  section,  prepares  the  satellite  data  for  height  deter¬ 
mination.  In  addition,  it  sorts,  screens,  and  combines  the 
conventional  data  received  from  the  Air  Force*s  global 
weather  data  base  to  produce  a  vertical  and  horizontal 
3DNEPH  gridded  layered  format. 

Next,  the  satellite  data  goes  to  the  visual  and 
infrared  processors.  The  VIS  processor  steps  are  outlined  in 
Fig.  14.  Each  pixel  grayshade  is  directly  compared  to  a 
grayshade  (background/surface  brightness)  specified  for  each 
25  n  mi.  square.  This  value  represents  a  clear  sky  earth 


radiance  for  each  square.  The  difference  between  each  pixel 


grayshade  and  the  background  brightness  is  used  in  Pig.  15 
to  decide  if  there  is  a  cloud  or  not.  The  curve  in  Pig.  15 
represents  the  brightness  threshold  for  clouds;  any  vari¬ 
ability  greater  than  this  is  a  cloud  area.  This  curve 
incorporates  the  background  brightness  and  the  varying 
brightness  associated  with  different  surfaces  (  e.g., 

desert).  The  total  number  of  cloud  decisions  divided  by  the 
total  number  of  available  pixels  represents  the  amount  of 
cloud  cover  in  the  square. 

The  infrared  processor  incorporates  a  slightly  dif¬ 
ferent  analysis  method.  The  major  steps  of  the  processor 
are  included  in  Fig.  16  as  a  flow  chart.  The  first  step  is 
to  construct  a  frequency  distribution  of  the  6«  raw  gray- 
shades  for  each  25  n  mi.  area  to  detect  significant  modes 
(maxima)  of  6  or  mere  grayshades  or  12  or  more  adjacent 
shades.  If  more  then  2  modes  are  identified,  the  modes  are 
combined  on  the  basis  of  relative  amplitude  (strongest  modes 
are  retained)  and  the  differences  between  modes.  If  there 
is  no  apparent  mode,  then  the  processor  defaults  to  a  cloud 
threshold  temperature  determined  from  appropriate  surface 
temperatures.  3efore  any  cloud  decision  is  made,  a  correc¬ 
tion  is  applied  to  the  IR  temperatures  to  compensate  for 
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atmospheric  absorption  and  other  effects  sach  as  zenith  and 
local  angle  which  create  biases  in  the  IS  data.  Pig.  17  is 
an  example  of  a  curve  used  to  derive  bias  correction  for  a 
particular  satellite.  The  actual  bias  correction  is  the  dif¬ 
ference  between  the  best  fit  curve  and  the  diagonal  line  of 
perfect  correlation.  Next  a  representative  temperature  is 
selected  for  each  mode  (usually  the  coldest  in  the  mode)  and 
these  data  are  compared  to  a  gridded,  conventionally  derived 
surface  temperature  data  base  to  sake  a  cloud/no  cloud  deci¬ 
sion.  Pig.  18  graphically  represents  this  decision  process. 
The  final  cloud  decision  is  made  by  using  Fig.  19.  This 
compensates  for  irregularities  in  the  surface  temperature 
data  base  which  limit  the  reliability  of  the  cloud/no  cloud 
decision.  The  cloud  amount  is  determined  by  the  number  of 
pixels  in  the  cloudy  mode  (Fye  and  Logan,  1977).  Addition¬ 
ally,  the  IH  processor  determines  the  cloud  height  by  com¬ 
paring  the  coldest  temperature  in  the  cloudy  mode  to  a 
real-time  upper  air  data  base. 

Other  processors  include  the  conventional  data  pro¬ 
cessor,  the  data  integration  processor,  and  the  manual  data 
processor.  The  conventional  data  processor  takes  the  sur¬ 
face,  aircraft,  and  upper  air  reports  and  produces  a 
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Piqure  17. 


An  IE  Bias  Correction  Curve 
Computer  Program  (from  Fye, 


horizontally  gridded,  eight  mesh  data  base  of  cloud  informa¬ 
tion  for  15  vertical  layers.  These  layers  of  data  and  the 
output  of  7IS  and  IE  processors  are  merged  by  the  data  inte¬ 
gration  processor.  The  processor  also  performs  meteorologi¬ 
cal  consistency  checks.  The  cloud  amounts,  heights,  and 
types  from  7IS  and  IB  data  are  converted  to  specific  cloud 
layer  amounts  and  heights.  Conflicts  are  also  resolved  by 
this  processor  by  checking  age  of  data,  greatest  cloud 
cover,  layer  amount,  and  cloud  cop  height.  Next,  the 


9 


NUMG 


Background  Brightness 
G  rayshade 


Figure  19.  ^il^d.CgtoJf^urv^Ose^jo  ^t#Low  clouds 


51 


satellite  and  conventional  data  are  compared  and  integrated 
with  a  continuity  field  to  complete  the  cloud  analysis 
(Fye,197  8).  The  continuity  field  is  the  previous  or  latest 
3DNEFH  analysis  which  fills  in  Hissing  data  and  parameters. 

The  analyst  next  enters  manually  encoded  input  to 
alter  the  3DNEPH  data  base  as  necessary.  The  manual  data 
processor  takes  this  information  and  incorporates  it  into 
the  overall  data  base.  This  provides  an  avenue  for  quality 
control,  correction  of  errors  and  improvement  of  the  overall 
analysis. 

3-  Biggegtral 

A  bispectral  method  for  determining  cloud  height  and 
amount  was  developed  by  Reynolds  and  Vonder  Haar  (1977). 
Visual  and  infrared  data  were  each  analyzed  quantitatively 
over  an  array  of  scan  spots.  This  analysis  yielded  cloud 
amount  and  cloud  top  temperature,  from  which  cloud  height 
was  derived  by  comparison  to  a  nearby  vertical  temperature 
sounding. 

The  method  was  based  on  manipulation  and  solution  of 
two  equations  (Eg.  1  and  2,  Table  IV)  .  These  equations  are 
basically  comprised  of  adjustments  for  the  fraction  of 
cloud-covered  and  clear-sky  areas  in  a  single  field  of  view 


and  for  backbody  emissivity  of  the  clouds 


One  was  the 


satellite  visual  radiance  equation  (Equation  1,  Table  IV) 

and  the  other  vas  the  satellite  infrared  radiance  equation 

(Eq.  2,  Table  IV).  The  equations  related  the  aaount  of  cloud 

(Acid  in  tenths)  to  either  cloud  albedo  (pcld  X  Hs)  in  the 

visual  or  cloud  radiances  (Held)  in  the  infrared  channels 

(note  that  Acid  plus  the  aaount  cf  clear  area,  Aclr  in 

tenths,  equaled  1.0).  The  IR  radiance  equation  also 

accounted  for  the  effect  of  emissivity  (e )  of  clouds  on 

cloud  and  surface  (Nclr)  radiances.  In  both  equations,  Aclr 

was  replaced  by  1.0-  Acid.  The  rest  of  the  terms  in  the 

equations  except  two  (Acid  and  Ncld)  were  available.  Average 

values  of  the  cloud  and  surface  albedoes  were  known  and  more 
# 

exact  values  can  be  obtained  in  practice  ."from  the  actual 
array  of  satellite  measurements”  (Reynolds  and  Vonder  Haar, 
1977).  Measured  surface  temperatures  were  converted  to  IR 
radiances  and  used  for  the  surface  radiance.  The  emissivity 
was  initially  assumed  to  be  0.9  for  all  clouds.  This  left 
only  the  two  unknowns  to  be  determined.  Eg.  1  (Table  IV)  was 
then  solved  for  Acid  and  Eq.  2  (Table  IV)  for  Ncld  radiance 
which  resulted  in  Eq.  3  and  4  (Table  IV).  The  Acid  value  was 
derived  first  from  the  visual  satellite  data.  Next,  Acid 
was  used  with  the  simultaneously  received  IR  data  to  obtain 


the  cloud  radiance  value.  The  resulting  radiance  was  then 
converted  to  cloud  top  teeperature  via  Plancks*  equation. 
The  final  step  was  to  use  a  local  upper  air  sounding  and  the 
cloud  temperature  to  obtain  the  cloud  top  height. 

Reynolds  and  Vonder  Haar  (1977)  performed  an  uncer¬ 
tainty  analysis  on  this  technique  to  dataraine  on  what  the 
uncertainty  (A)  of  the  cloud  amount  (Acid)  and  cloud 
temperature  (dNcld)  measurements  depend  and  how  large  an 
uncertainty  (A)  is  allowable.  It  was  shown  that  A  icld 
depends  on  the  m  agnitude  and  precision  of  Hs,  Ms,  pcld,  and 
pclr,  where  A  Hs  and  A  pc  lr  will  be  less  than  10*  of  the 
error.  Ncld  depends  on  the  magnitude  and  uncertainty  of  Ml, 
Acid,  Nclr  and  e. 

The  two  quantities  A  Acid  and  dNcld  were  evaluated 
using  typical  mid-latitude  values  obtained  from  satellites. 
Table  V  gives  the  absolute  dAcld  error  for  given  A  Hs,  Ms, 
and  pcld.  These  results  show  that  the  "uncertainties  must 
stay  below  the  10 %  level  if  reasonable  accuracy  is  to  be 
maintained”  (Reynolds  and  Vonder  Haar,  1977)  .  This  accu¬ 
racy  is  possible  through  use  of  ground  and  cloud  truth  meas¬ 
urements  along  with  improvements  in  instrument  calibration. 
At  present  there  is  approximately  a  3  degree  Celcius  (°G) 
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Algorithas  for  Cloud  Top  Temperature  and  Cloud  Amount  (afte 
Reynolds  ana  fonder  Haar,  13  77) 


(1)  *s 
(2) 

<3>  Acld 


AclrpclrHs  *  Acld#c1dHs 

Ac1r  Nc1r  *  Acld  Nc1d  * 

*  Acld(l*°  ’  «>"clr 

ms  -  9c1r  Hs 
Hs^scld*  ®clr^ 


<«>  "eld 


"L  -  Nclr 
*  Ac1d 


t-  N 


clr 


(5) 

(6) 


Acld  +  Ac1 r 


c 


1.00 

Nc1r 


N 


dr 


cld  •  cloud 
dr  ■  cl  ear 

A„, j  »  area  of  array  covered  with 
c1d  cloud 

«s  *  measured  visible  radiance 
Bcld  *  b1_<1'*rect*onal  reflectance 
cloud  areas 

H.  *  solar  irradiartce  reaching 
"*  the  surface  (  .5  <  x  < 

.8  um) 

p  .  -  bi-directional  reflecting 

c  r  of  clear  areas 

Nr1(4  *  spectral  !R  (10  <  x  < 

C1Q  12  um) 

N  1  «  spectral  IR  radiance  from 

clear  area 

-  measured  IR  radiance 
c  *  IR  emissivity 
A.,-  •  area  of  array  clear  of 


error  in  satellite  derived  teeperature  which  is  equivalent 
to  approxisately  5%  A  in  Nclr.  Held  will  have  a  16°C  er¬ 
ror  if  the  Acid  has  24%  error  (fable  V) .  This  error  can  be 
reduced  by  better  Acid  detersina tion  using  higher  resolution 
data  and  through  use  cf  ground  truth  seasuresents. 

TABLE  V 

The  Errors  in  Cloud  Aaount  Jwith  dHs  *  5%)  and  in  Cloud 
Teeperature  (Held  =  265  Kr  after  Heynolds  and  vonder  Haar. 

1977) 


Uncertainty 


A  Ac  Id  (Acid 
AAcld  (Acid 
A Ac Id  (Acid 


Ncld 
AZ 


5Z  p clr, 
pcld,  Ms 


IlHjI 

n 


10Z  pclr, 
pcld  Ms 


15Z  pclr 
Pcld,  Ms 


i  0.14 

0.20 

0.21 

0.30 

0.29 

0.40 

5Z  ML, 

10Z  ML, 

Nclr,  e 

Nclr,  e 

24Z  in  Acid 

40Z  in  Acid 

30° 

4.2  km 


The  bispectral  method  was  tested  with  NOAA  satellite 
coverage  in  an  area  that  had  cloud  observations  and  upper 
air  data  available.  Three  locations  were  chosen;  White  Sands 


Missile  Sange  (WSMB)  ,  Denver,  and  Oklahoma  City 


A  75  X  75 


km  area  was  used  because  it  best  corresponded  to  surface  and 
upper  air  data.  Average  values  of  pcld,  pclr,  and  Hs  were 
assuaed  initially  and  then  changed  as  sun  angle  and 
satellite  viewing  angles  varied  (Reynolds  and  Ponder  Haar, 
1977).  The  results  showed  a  res  error  of  0.2  in  cloud 
amount  with  a  slight  underestimation  bias.  The  cloud  height 
error  with  cirrus  included  was  4.6  km  with  a  -4.2  km  bia3. 
without  cirrus  cloud,  the  rms  error  improved  to  0.5  km  with 
overestimation  of  0.27  km. 

This  problem  with  cirrus  is  linked  to  the  use  of  0.9 
for  the  emissivity  of  all  plouds.  Because  of  the  ice  crys¬ 
tals  in  cirrus,  its  emissivity  is  less  than  the  0.9  value. 
Two  approaches  were  considered  by  Reynolds  and  Vonder  Haar 
(1977)  to  correct  this  problem. 

One  approach  was  to  correlate  the  visible  spectral 
albedo  of  a  cloud  to  its  emissivity  through  the  blackbody 
temperature.  Pig.  20  shows  the  scatter  diagram  of  the  rela¬ 
tionship  between  the  albedo  and  the  blackbody  temperature 
and  a  second  order  fit  which  was  used  to  determine  the  cloud 
top  temperature  for  a  specific  albedo.  Once  this  cloud  top 
temperature  is  known  the  emissivity  can  be  calculated 
using  Eg.  4,  Table  17  (derived  from  Eg.  2  of  Table  17). 
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Fig.  21  is  the  resulting  curve  which  relates  albedo  to  eais- 
sivty.  The  curve  is  a  second  order  equation  of  the  best  fit 
curve  seasured  in  the  bispectral  sethod  along  with  a  slight 
adjustsent  to  pclr  and  pcld.  This  cut  the  error  of  the  old 
aethod  in  half  (Reynolds  and  yonder  Baar,  1977). 

The  second  procedure  for  deteraining  cirrus  cloud 
height  was  to  coapare  channels  4  and  8  of  the  Vertical 
Teaperature  Profile  Radioaeter  (VTPR)  on  the  NOAA-4 
satellite.  This  led  tc  a  reduction  by  one  half  in  the  appar¬ 
ent  error  of  the  old  aethod.  Table  VI  lists  the  results  of 
using  a  cirrus  correction  and  the  VTPR  approach  compared 
with  non-ad justaent  for  cirrus. 


D.  PRECIPITATION 

Four  approaches  for  esti sating  precipitation  from 
satellites  are  reviewed  in  this  section.  They  are  cloud 
index  (Barrett  and  Martin,  1981;  Barrett,  1981),  life  his¬ 
tory  (Barrett  and  Martin,  1981;  Griffith  and  Woodley,  1981), 
bi-spectral  (Liljas,  1981a,  1981b;  Austin  and  Lovejoy, 

198  1),  and  cloud  modeling  (Wylia  ,  1979;  Simpson  and  Wiggert, 
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Fiqure  21.  The  Relationship  Between  Albedo  and  IR 

Emissivity  from  Empirical  Studies  (from  Reynolds 
and  Vonder  Haar,  1977). 


TABLE  VI 

Comparison  of  cirrus  correction  and  VTPR  results  (from 
Reynolds  and  Vonder  Haar,  1977) 


CIRRUS  COWECriM  RESULTS 


TTN  RESULTS 


R6  Error  •  1.7S  ka 
lltl  ■♦l.i  ka 


i.  sigaa  mass 

Barrett  (1981)  developed  a  cloud  index  method  of 
deteraining  rainfall  (called  the  'Bristol  Method').  It  is  a 
manual  aethod  which  uses  observed  rain  gage  data  and  the 
cloud  index  system  (see  Section  B)  to  pradict  the  aaount  of 
rainfall  accumulated  over  a  specified  tiae  period  (usually 
12  hours)  in  grid  squares  of  a  selected  size  (1/6°  to  1° 
grid  squares).  Cloud  area  and  altitude  are  also  taken  into 
account  (Barrett  and  Martin,  19  81).  Fig.  22  presents  the 
flow  diagram  of  a  recent  version  of  the  Bristol  Method  with 
two  stages,  preparatory  and  operational.  Both  IR  and  VIS 
are  used  whenever  possible  to  identify  cloud  types.  Hext, 
the  cloud  types  of  each  grid  sguare  are  assigned  an  index 
(see  Table  VII).  The  cloud  index  is  generally  "evaluated 
differently  for  ranges  of  chosen  cloud  types  and  for  differ¬ 
ent  regions"  (Barrett  and  Martin,  1981).  Then  mean  cloud 
indices  are  calculated  along  with  mean  observed  rainfall 
amounts  for  each  cloud  field  unit.  Subsequently,  regression 
techniques  are  used  to  relate  the  observed  rainfall  to  the 
cloud  type,  and  these  help  establish  rainfall  estimates  for 
areas  of  satellite  cells  where  conventional  data  are  sparse 
(see  Pig.  23).  Tests  show  that  the  accuracy  of  this  method 
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Figure  22.  Flow  Chart  for  a  Recent  Version  of  the  Manual 
Nephanalv sis  (from  Barrett  and  Martin,  1981). 


TABLE  7 II 


Empirically, derived  Rainfall  Probabilities  and  Intensities 
from  Satellite  Pictures  (  from  Barrett  and  Martin,  1981) 

1  2  3 

Assigned  probabili-  Assigned  intensities 

States  of  the  sky  ties  of  rainfall  of  rainfall  (rela- 


(nephanalysis  cloud 
categories) 

(relative  scale  range 

0  -  1.00) 

tlve  scale  range 
0  -  1.00) 

Cumulonimbus 

0.90 

0.80 

Stratiform 

0.50 

0.50 

Cumulif orm 

0.10 

0.20 

Stratocumulif orm 

0.10 

0.01 

Cirriform 

0.10 

0.01 

Clear  skies 

«... 

2.  Life  History 

The  Grif fith/Woodley  (Griffith  and  Woodley,  1981; 
Griffith  et  al.,  1973)  technique  for  estimating  convective 
precipitation  was  developed  using  a  cloud  life  history  model 
in  which  satellite  imagery  is  compared  to  radar  echo  life 
histcry  of  the  cloud  area  (Griffith  and  Woodley,  1981;  Bar¬ 
rett  and  Martin,  1981).  The  basis  of  the  technique  is  the 
hypothesis  that  in  the  tropics  rain  occurs  with  cold 
(bright)  clouds  and  the  intensity  of  the  rain  is  a  function 
of  the  stage  of  development  of  cumulonimbus  clouds.  In  a 
study  of  clouds  over  Florida,  Griffith  and  Woodley  (1981) 
found  that  the  cumulonimbus  cloud  area  is  related  to  radar 
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Figure  23. 
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echo  area  (see  Fig.  24).  Cloud  areas  (Ac)  are  evaluated  in 
either  Vis  or  IR  satellite  imagery  and  then  compared  to  the 
radar  echo  area  (Ae)  .  Fig.  24  is  a  plot  of  the  relationship 
between  cloud  and  echo  areas  (Griffith  at  a^. ,  1978).  it  was 
empirically  derived  using  measurements  and  comparisons  of 
Florida  clouds  viewed  by  radar  and  geosynchronous 
satellites.  The  rain  cloud  threshold  (IF)  used  in  the 


approach  was  25  3  K 
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Figure  24. 


Cloud 


normalized  to  relative  maximum  cloud  area  (Am) 
(from  Griffith  £t  a|,.,  1978). 


This  technique  is  fully  automated  (Figure  25)  and 
has  been  tested  in  the  Florida  area  and  in  extra-tropical 
convective  areas  (Barrett  and  Hartin,  1981).  There  are  four 


major  computer  modules  that  together  produce  a  rain  map.  The 
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programs  are  navigation,  cloud  isolation  and  tracking,  cloud 
life  history  compilation  and  volumetric  calculation  of  rain, 
and  rain  napping  (Griffith  and  loodley,  1981;  Barrett  and 
Martin,  1981)  • 

This  nodal  has  been  applied  on  several  tine  scales 
and  in  different  regions  as  shown  in  Table  7III.  The  error 
statistics  (bias  (B)  ,  root  mean  square  error  (Bras)  ,  and 
linear  least  squares  fit  (R,  slope,  and  intercept))  for 
these  applications  (Table  IX)  show  that  on  average  the  nodel 
results  agree  fairly  well  with  the  ground  truth  data 
(Griffith  and  Woodley,  1981).  The  Eras  is  small  except  for 
two  locations,  with  an  overall  average  of  9.35  millimeters 
per  hour.  The  bias  is  close  to  1.0  0  which  is  the  value  of 
perfect  correspondence.  In  the  linear  fits,  all  correla¬ 
tions  (R)  are  greater  than  or  equal  to  9.78,  where  1.0  is 
perfect  correlation. 

Lovejoy  and  Austin  (1979)  developed  a  technique  in 
which  visual  and  infrared  wavelength  patterns  are  used  to 
identify  clouds  with  a  high  probability  of  rain.  These 
clouds  have  large  optical  and  physical  thickness  "as 
measured  by  the  visible  albedo  in  comparison  with  their 
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Figure  25.  The  Griff  ith/Wcodley  Technigua  Flow 
for  GATE  (from  Barrett  ana  Sartrn. 
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TABLE  ?III 


Appncations^t^GriJfJthgooU^^ 


;v.uua^uic  v  —  —  w  ' 


AREA 

Global  Atmospheric  Research  Project' * 
Atlantic  Tropical  Experiment  (GATE) 

Florida  Area  Cumulus  Ex per inert 
TFACE) 

equatorial  Pacific  Ocean  Climate 
Studies  (EPOCS) 

High  Plains  Experiment  (HIPLEX) 

Hurricanes 
Flash  floods 


PERIOO 

June  27  to  Sept.  20,  1974 


Selected  days  June  to  September 
1972  to  1980 

2-1/4  days  November ,  1978 
30  days  November,  1979 

Selected  days.  May  to  July 
1976  to  1978 

Selected  storms  1969  to  present 

Biq  Thompson,  CO  7/31  -  8/1,  1976 
Johnstom,  PA  July  19-20,  1977 


tl 
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TABLE  XX 


The  Error 


tatistics  fro®  the  Applications  List? 
Ill  (from  Griffith  and  Woodley,  1981) 


Listed  in  Table' 


LOCATION 

N 

B 

Sms 

R 

SLOPE 

INTERCEPT 

GATE 

53 

8.84 

0.79 

0.87 

0.97 

-1.47 

HIPLEX 

(dense) 

15 

0.77 

2.81 

8.90 

0.62 

0.07 

FACE 

11 

0.80 

0.42 

0.95 

1.18 

-0.23 
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height,  determined  by  the  intensity  of  the  IB  emmission" 
(Austin  and  Lovejoy,  1981).  Liljas  (1981a,  1981b)  followed 
this  with  a  simliar  type  of  bispectral  identification  of 
precipitation  area  based  on  7IS  and  IR  iiagery.  This  sub¬ 
section  will  first  review  Austin  and  Lovejoy's  (1981) 
technique  and  then  Liljas1-  (1981a,  1981b)  technique. 

Austin  and  Lovejoy  (1979)  first  addressed  the  ques¬ 
tion,  "How  well  can  rain  areas  be  delineated  in  coincident 
visible  and  infared  satellite  images?"  They  felt  that  the 
problem  was  to  find  boundaries  between  rain  and  no  rain. 
Estimation  of  rain  amount  is  assumed  to  depend  on  two  steps, 
the  delineation  of  the  rain  areas  and  estimation  of  rain 
rate  once  the  area  is  known.  The  authors  used  three  statis¬ 
tics  to  characterize  errors  and  from  these  were  able  to 
estimate  rain  area  and  rain  rate.  The  first  step  in  rain 
estimation  was  to  map  a  frequency  matrix  for  both  rain  areas 
and  non-rain  areas  (see  Fig.  26  and  27)  .  The  ratio  of  rain 
to  non-rain  was  found  and  plotted  (see  Fig.  28)  along  with  a 
line  of  greater  then  50%  probability  of  rain.  From  this,  a 
rain  map  (see  Fig.  29)  was  produced.  The  next  step  was  tc 
try  to  quantitatively  estimate  the  rain  rate.  Austin  and 
Lovejoy  (1981)  tried  several  techniques  to  do  this  but  met 
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with  success.  They  believe  this  was  due  to  "...the 

visible  and  IS  wavelengths  primarily  responding  to  the  rela¬ 
tive  abundance  of  cloud  droplets  and  not  to  precipitation 
sized  particles"  (Austin  and  Love  joy,  1981)  .  They  concluded 
that  a  good  or  even  passable  guantitativa  estimate  of  rain¬ 
fall  rate  is  not  realistic  or  practical  using  visual  and 
infrared  data.  The  authors  felt  their  study  supported  the 
hypothesis  that  visual  and  infrared  satellite  images  are 
good  for  estimating  rain  areas  but  not  for  estimating 
quantitative  rain  rates. 

Liljas*  (198  1a)  main  objective  was  to  apply  an 
approach  similiar  to  that  used  by  Lovejoy  and  Austin  (1979) 
for  identifying  precipitation  areas  and  intensity  using 
TISOS-S  satellite  imagery,  but  without  using  radar  to  cali¬ 
brate  it.  He  decided  not  to  try  quantifying  intensity,  but 
to  instead  define  it  qualitatively  by  categories  of  light  to 
heavy  rain.  This  fellows  the  work  of  Muench  and  Keegan 
(1979)  on  the  relation  of  cloud  reflectance  and  intensity  of 
rain.  Pursuing  their  concepts,  Liljas  (1981a,  1981b)  devel¬ 
oped  a  relative  precipitation  intensity  identification 
technique  which  relies  on  his  cloud  typing  model  (see  Sec¬ 
tion  B,  Chapter  II) .  It  can  be  applied  quickly  no  both 
frontal  and  convective  type  precipitation  and  gives  the 
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Figure  26. 


Frequency  Distributi 
for^the  No-rain  Casa 


on  of  ?is  and  IB  Data  Pal 
(from  Barrett  and  Martin 


Figure  27 


IVjtt  HL  •  W  «  «/ 
i  i!Ny  m  n  it  Mtv 


•Ml  •  • 

•  • 
Nt  VJ  i 


Brightness  -* 


Figure  28.  Probability  of  fiain  i$  Percent  from  the  Arrays 
of  Fig.  2o  and  2/  with  the  50%  Boundary 
Sketched  (froa  Barra  tt  and  Martin,  1981). 
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Figure  29. 


The  Resulting  Satellite  Sain  Sap  froa  the  IJse  o 
♦  ne  Saae  Cats  as  Fig.  26  ana  27.  The  vertical 
lines  are  radar  pain  areas  and  horizontal  lines 
are  satellite  ram  areas  (hatched  areas  are  the 
agreement  areas)  (from  Barrett  and  Martin, 
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information  in  a  real-time  fashion.  The  first  step  is  to 
identify  the  cload  types  (see  Fig.  9).  Best,  the  aiea  on  the 
graph  representing  ccaulonimbus  and  nimbostratus  is  sec¬ 
tioned  (see  Fig.  10)  based  on  the  preeise  that  cold,  high 
and  dense  clouds  give  the  highest  precipitation  rate.  in 
additional  requirement  for  the  cloud  top  temperature  to  be 
colder  than  -22°C  was  established  (Liljas,  1981a)  and  sup¬ 
ported  by  invest igaticns  of  Muench  and  Keegan  (1979)  .  Table 
X  gives  the  thresholds  determined  by  summing  the  infrared 
and  visual  digital  values  that  correspond  to  the  precipita¬ 
tion  intensity  levels.  The  method  gave  good  information  on 
the  distribution  of  precipitation  areas  and  the  relative 
intensity  within  the  areas  for  both  frontal  and  scattered 
showers.  In  addition,  the  method  detected  "intense  precipi¬ 
tation  in  mesoscale  cloud  clusters  between  synoptic" 
observations  (Liljas ,  198  1  a) . 

4.  Cloud  Model 

Wylie  (1979)  compared  tropical  and  mid-latitude 
cloud  area  and  rain  rates  to  see  if  the  Griffith/Woodley 
(Griffith  et  al«#  1978)  and  other  techniques  could  be 
applied  to  non-tropical  cloud  features,  instead  of  the  trop¬ 
ical  cloud  systems  for  which  they  were  developed.  Six  days 


TABLE  I 


Chi  +  Ch4  291-310  light  rain 

311-330 
331-350 
351-370 
371-390 

>390  very  strong  rain 


of  t mo  different  types  of  rain  clouds  (frontal  convective 

and  frontal  stratiform)  in  the  Montreal,  Quebec,  Canada  area 

were  chosen  for  comparison  to  the  1979  GATE  data  set  of 

0 

Stout  et  ajj,.  (1979)  (see  Table  71 1)  . 

The  rain  rates  of  the  two  areas  are  shown  in  Fig. 
30.  The  GATE  rain  rates  are  higher  than  the  rain  rates  of 
Montreal,  even  at  the  colder  thresholds.  The  Montreal  rain 
rates  tend  to  vary  considerably,  especially  the  convective 
cases  (June  cases)  .  This  shows  that  a  scheme  based  on  the 
GATE  data  would  overestimate  rain  rate  in  the  midlatitudes, 
and  additionally,  if  one  uses  a  constant  adjustment,  there 
would  be  a  large  variance  in  the  rain  amounts.  Table  XI 


shows  the  comparison  for  GATE  and  Montreal  rates  of  precipi¬ 
tation.  Wylie  (1979)  found  that  the  precipitable  water  -r- 
alleled  the  difference  in  the  average  rain  rate  with  sose 
minor  exceptions  when  the  troposphere  was  exceptionally  sta¬ 
ble.  He  decided  to  use  Simpson  and  Wiggert's  (1969)  one-di¬ 
mensional  cloud  model  to  assess  the  effects  of  stability  and 
to  possibly  improve  his  results  in  the  mil-latitudes. 


Intarvd  can***  M  (K) 

Figure  30.  Rain  Rates  Averaged  Over  Cloud  Areas  Measured  on 
IR  Images  for  Various  £loud  Threshold 
Temperatures  (from  Wylie,  1979). 

This  model  simulates  cloud  growth  by  using  an 
ascending  2  km  bubble.  The  model  incorporates  cloud  physics 
processes  in  the  development  of  precipitation  (Simpson  and 
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wiggert,  1969)  .  Only  one  bubble  ascent  was  used  although 
normally  the  babble  would  ascend  nany  tines.  &s  a  result, 
the  rainfall  prediction  could  not  be  used  directly,  but  was 
nevertheless  used  as  an  adjustment  to  correct  the  Montreal 
data  for  stability. 

Table  XI  shows  the  results  of  Wylie's  (1979)  study: 
the  corrected  precipitation  of  the  satellite  closely  matches 
the  radar  neasured  amount  of  rain.  He  noted  that  the  over* 
all  bias  was  small,  and  therefore  by  using  a  stability  cor* 
rection  (calculated  from  the  on e-dimensional  cloud  model), 
the  technique  to  measure  precipitation  by  satellite  in  the 
tropics  could  be  applied  to  mid- latitude  areas. 


hi.  sim  ciggg  gam  gfissaiEiifii 

1.  INTRODUCTION 

The  cloud  analysis  aodel  adapted  for  the  SP1DS  interac¬ 
tive  coaputer  is  a  combination  of  procedures  drawn  froa  sev¬ 
eral  previously  tested  techniques.  The  mainstay  of  the 
model  is  Liljas'  (19  81a)  paral lelpipedon  (Fig.  9  and  10) 
classification.  added  tc  this  are  a  texture  test  for  dis¬ 
criminating  some  cloud  types,  a  non-linear  least  squares 
curve  fit  to  discern  cirrus  and  alto  stratus,  Harris  and  Bar¬ 
rett's  (1978)  cloud  amount  calculation,  part  of  Reynolds  and 
vonder  Haar's  (1977)  cloud  top  temperature  calculation,  and 
Liljas'  (1981a)  qualitative  precipitation  intensity  model. 
The  following  sections  describe  these  parts  of  the  SPADS 
model  and  explain  why  they  were  selected.  The  last  section 
of  the  chapter  summarizes  with  the  flow  chart  of  the  model's 
computer  program. 

B.  CLOUD  TYPING 

The  approach  selected  for  cloud  typing  was  developed  and 
tested  by  Liljas  (1981a).  It  is  a  aultispectral  method 
that  uses  three  channels  to  type  clouds  and  differentiate 
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between  land  and  water 


The  SPADS  cloud  aodel  uses  two 


channels,  thereby  becoming  bispectral  (the  near-IB  channel 
was  not  used  due  to  non-availability  in  GOES  VISSR  data). 
Fig.  10  depicts  the  bispectral  aodel  incorporated  in  the 
SPADS  Cloud  Model.  Both  the  visual  and  infrared  thresholds 
shown  in  Fig.  10  were  converted  froe  Tiros-N  Advanced  Very 
High  Resolution  Radiometer  (A  VRRB)  to  GOES  VISSR  data 
counts.  Liljas*  visual  data  brightness  counts  ranged  from  0 
to  255  while  the  SPAOS  counts  range  from  3  to  63.  A  four  to 
one  linear  relationship  exists  in  the  VIS  channel  and  there¬ 
fore  a  direct  conversion  of  thresholds  was  performed.  The 
IR  thresholds  were  converted  to  temperatures  (K)  using  an 
August  1980  TIROS-N  AVHRR  tea perature  calibration  table. 
This  gave  a  rain  cloud  threshold  of  about  251  K  which  is  in 
agreement  with  the  Griff ith/Woo dley  (Griffith  and  Woodley, 
1981)  and  the  Liljas  (1981a)  rain  cloud  thresholds.  These 
temperature  thresholds  were  converted  to  GOES  IR  brightness 
counts.  The  transformed  thresholds  used  in  the  SPADS  cloud 
model  are  depicted  in  Fig.  31. 
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One  of  the  reasons  for  selecting  the  Liljas  eodel  is 
that  it  was  tested  or  an  interactive  computer  siailiar  to 
SPADS.  In  addition*  the  test  area  selected  by  Liljas 
(1981a)  included  land  and  ocean  areas  along  with  nany  types 
of  clouds  associated  with  oesoscale  and  synoptic  scale  fea- 
tures.  A  complete  description  with  color  pictures  is  in 
Liljas' (  1981a)  technical  report.  In  addition,  a  combined 
approach  using  texture  and  bispectral  counts  was  suggested 
by  Liljas  (1981a),  Harris  and  Barrett  (1978)  and  Reynolds 
and  Vonder  Haar  (1977).  The  U.S.  Air  Force's  3DNEPH  program 
used  both  these  approaches,  but  documentation  of  thresholds 
and  tests  was  not  available  to  us.  Thus,  the  3DNEPH  program 
could  only  be  used  as  a  qualitative  guide. 

The  GOES  IR  resolution  is  approximately  4  n  mi.  at  the 
subsatellite  point,  but  for  the  visual  channel  a  higher  res¬ 
olution  of  0.5  n  mi.  is  available.  This  Vis  resolution 
affords  an  opportunity  to  get.  a  better  estimate  of  cloud 
cover.  All  the  studies  reviewed  highly  recommended  the  use 
of  an  array  of  brightness  counts  to  give  an  improved  repre¬ 
sentation  of  the  cloud  picture.  The  0.5  r.  mi.  resolution 
allowed  an  8X8  array  of  visual  pixels  per  IS  pixel.  For 
each  grid  space  an  average  visual  brightness  count  is 
calculated  (Eg.  1  ,  Table  II). 


At  the  saae  time  as  the  visual  average  brightness  is 
being  calculated,  the  standard  deviation  of  the  cloudy  vis¬ 
ual  counts  per  grid  scace  is  determined.  The  standard  devi¬ 
ation  (Eg.  4,  Table  I)  is  the  statistic  used  to  represent 
texture.  As  suggested  and  tested  by  Harris  and  Barrett 
(1978),  the  standard  deviation  is  used  to  discriminate 
between  stratiform  and  cumulus  humulus,  small  and  large 
cumulus  congestus,  and  altostratus  and  cirrus.  The  initial 
standard  deviation  threshold  values  for  this  discriminant 

analysis  were  approximated  by  using  information  in  Fig.  6 
0 

and  8.  An  additional  test  was  used  first  for  the  cirrus/al- 
tostratus  decision.  A  second  order  polynomial  least  squares 
fit  was  applied  to  the  cirrus  (c)  points  and  also  to  the 
altostratus  (d)  points  in  Fig.  10.  The  calculated  curves 
are  depicted  in  Fig.  31,  with  the  associated  equations  in 
Table  XII.  Entering  the  equations  with  the  average  visual 
brightness,  the  resulting  13  value  is  compared  to  the  actual 
13  to  determine  cloud  type.  If  the  calculated  value  falls 
between  the  two  curves  then  an  additional  test,  using  the 
standard  deviation  threshold,  is  used  to  discriminate  cloud 
type. 
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C.  CLOUD  TOP  TEMPER* TUBBS 

The  Reynolds  and  Yonder  Haar  (1977)  study  is  used  as  a 
basis  for  estimating  cloud  top  teaperatures.  Their 
calculation  is  based  cn  two  equations  (Eq.  1  and  2,  Table 
IT) ,  one  for  cloud  aaount  and  the  other  for  the  cloud  top 
temperature.  The  resolution  of  the  data  used  was  4  n  ai. 
but  the  authors'  analysis  suggested  an  improved  resolution 
would  give  better  results.  As  nentioned  above,  this  led  to 
the  choice  of  0.5  n  a  i.  visual  data  in  this  study. 

The  first  step  in  determining  the  cloud  top  temperature 
is  to  calculate  the  aaount  of  cloud  per  IR  pixel.  The  cloud 
decision  was  based  on  comparing  the  visual  digital  count  of 
each  pixel  in  the  8X8  grid  to  a  no  cloud  threshold.  This 
threshold  value  (initially  5  VI S  counts)  was  derived  from 
Fig.  31  (after  Liljas,  1981a)  where  it  is  the  lowest  VIS 
count  to  have  clouds.  A  cloud  threshold  test  was  also  used 
by  Harris  and  Barrett  (1978).  Their  threshold  value  was 
derived  by  comparing  the  visual  image  with  a  computer  pro¬ 
duced  density  slice,  but  neither  the  value  nor  the  formula 
to  obtain  the  value  were  included  in  their  report.  The  U.S. 
Air  Force's  3DMEPH  program  based  the  cloud  decision  cn  the 
aaount  of  visual  grayshade  difference  between  a  3  n  mi. 


pixel  and  a  background  brightness  (Fig.  15,  refer  to  Chap- 

« 

ter  II  for  nore  cob p let e  details).  This  prograa  uses  an 
extensive  global  background  brightness  file  that  is  updated 
regularly  (Fye,  1978).  The  developaent  and  saintenance  of 
such  a  file  for  SP&DS  is  not  feasible  at  this  tine. 

The  average  cloud  anount  is  determined  by  summing  the 
number  of  cloud  decisions  in  the  grid  space  and  dividing  by 
the  total  number  of  pixels  per  grid  space.  This  calculation 
followed  the  Harris  and  Barrett  (1  978)  Sq.  3  (Table  I). 
This  cloud  amount  fraction  is  then  used  in  the  Reynolds  and 
Vonder  Haar  (1977)  equation  for  cloud  top  IR.  radiance  (  2q. 
2,  Table  IV).  On?  of  the  reasons  this  cloud  amount 
calculation  is  used  instead  of  Reynolds  and  Vonder  Haar's 
equation  was  that,  in  the  latter,  three  additional  constants 
would  have  to  be  calculated  and  verified,  which  is  beyond 
the  scope  of  this  thesis.  The  authors  felt  that  it  was 
essential  in  using  their  cloud  amount  equation  that  these 
constants  be  well  determined  (Reynolds  and  Vonder  Haar, 
1977)  . 

The  next  parameter  of  the  cloud  top  radiance  equation  to 
be  determined  is  the  clear  area  IR  spectral  radiance.  Fleet 
Numerical  Oceanography  Center  ( FNOC)  analyses  of  surface 


temperatures  (K)  are  obtained  and  a  representative 
teaperature  profile  witching  the  location  of  the  IH  isage 
center  point  is  selected.  These  temperatures  are  converted 
to  IB  radiances  via  Planck's  function  and  used  -in  the  cloud 
top  teaperature  equation. 

Eaissivity  of  the  clouds  also  oust  be  deter  mined. 
Reynolds  and  Vonder  Haar  (1977),  among  others,  found  that 
the  eaissivity  is  nearly  constant  (approximately  equal  to 
0.9)  for  all  clouds  except  cirrus,  for  which  it  is  approxi¬ 
mately  0.55.  In  the  program,  cirrus  clouds  are  identified  by 
the  cloud  typing  and  standard  deviation  lodules  so  that  the 
lower  eaissivity  can  be  used  in  the  cloud  top  teaperature 
equation.  In  addition,  an  eaissivity  of  1.0  was  used  for 
cumulonimbus  and  niafcostratus  based  on  the  assumption  that 
these  clouds  do  not  allow  any  surface  radiance  through. 

The  cloud  top  radiance  (Ncld)  is  calculated  using 
Reynolds  and  Vonder  Haar's  (1977)  equation  (2q.  2,  Table 

III)  which  consists  of  the  surface  temperature  field 
obtained  from  FNOC  (converted  to  an  IR  radiance  via  Planck's 
law) ,  the  emissivity  (cirrus  or  non-cirrus) ,  the  satellite 
IR  radiance  (HI)  value,  and  the  fractional  amount  of  cloud 
(&cld)  .  Once  Ncld  is  calculated,  a  conversion  to  a  cloud  top 
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temperature  is  made  via  the  inverse  of  Planck's  function. 
This  teaperature  can  then  be  coapared  to  the  latest  PHOC 
upper  air  data  analyses  to  find  the  corresponding  height  in 
aillibars  (ab)  of  the  cloud  top.  The  upper  air  data  are 
available  at  ten  standard  levels  consisting  of  1000,  900, 

850,  700,  500,  400,  300,  250,  150,  and  103  mb. 

D.  PRECIPITATION  INTENSITY 

The  precipitation  intensity  identification  technique 
selected  is  the  Liljas*  (1981a)  method.  It  is  an  extension 
of  his  cloud  typing  method  (sea  Chapter  II)  and  therefore 
easily  fits  into  the  overall  SPADS  cloud  model  without  add¬ 
ing  much  time  to  the  process.  Wylie's  (1979)  cloud  model,  as 
discussed  in  Chapter  H,  used  the  Griffith/Woodley  (Griffith 
et  sJ,.,  1978)  technique  with  a  stability  correction  which 

gave  good  indication  cf  rain  rates.  One  problem  with  it  is 
that  the  Griffith/Woodley  technique  is  based  on  the  develop¬ 
ment  and  identification  cf  convective  cloud  systems'  life 
history  through  use  cf  radar  echoes  not  available  to  this 
study.  Another  problem  with  using  this  scheme  is  the 
computer  run  time  it  would  add  to  the  SPADS  cloud  model. 

Although  time  is  an  important  factor,  it  is  not  the  only 


reason  for  selecting  the  Liljas  method, 
reported  good  success  with  thi3  technique. 


Liljas  (1981a) 
obtaining  good 


indications  of  relative  precipitation  over  land  and  sea 
areas  as  verified  by  synoptic  reports.  He  felt  that  this 
technique  would  be  useful  for  sea  areas  where  no  synoptic 
observations  were  available.  In  addition,  relative  precipi¬ 
tation  from  aesoscale  cloud  clusters  was  detected  between 
synoptic  observations.  This  identification  scheme  was  modi¬ 
fied  for  use  in  the  SPADS  cloud  model  by  reducing  the  number 
of  precipitation  intensity  categories  fro*  six  (see  Pig.  9) 
to  three.  Fig.  32  depicts  the  new  categories  of  precipita¬ 
tion  types  representing  light,  moderate,  and  heavy  rainfall 
and  Table  XIII  includes  the  equation  for  the  lines  separat¬ 
ing  the  categories.  The  light  rainfall  category  includes 
the  requirement  that  cloud  top  temperatures  be  less  than 
251  K. 

E.  FLOW  CHAHT 

The  SPADS  Cloud  Hodel  (Fig.  33)  begins  with  obtaining 
the  initial  infrared  and  visual  data  fields  from  the  GOES 
satellite  and  conventional  surface  and  upper  air  data  from 
FNOC.  The  IS  and  VIS  satellite  data  are  collocated  both  in 
space  and  time.  The  upper  air  and  surface  temperatures  are 
selected  to  correspond  to  the  I,J  position  closest  to  the 


center  point  of  the  image  from  the  latest  FNOC  analysis  on 
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the  sane  date  as  the  satellite  pictures.  Next,  the  SPADS 
cloud  model  computer  program  (see  Appendix  A)  is  run. 

As  shown  in  Fig.  33,  the  first  step  of  the  prograa  is  to 
calculate  three  values  based  on  the  8X8  grid  of  VIS  pixels; 
the  average  VIS  brightness,  standard  deviation,  and  cloud 
amount  as  discussed  in  Section  3  of  this  chapter.  Once  the 
values  are  obtained,  the  cloud  typing  module  (based  on 
Liljas,  1981a)  is  called.  Two  tests  are  performed  within 
this  module;  one  is  based  on  comparison  of  the  IB  and  aver¬ 
age  Vis  counts  and  the  other  supplements  the  first  by  using 
texture  (standard  deviation)  to  complete  the  cloud  typing 

task.  If  nimbostratus  or  cumulonimbus  clouds  are  identi- 

# 

fled;  the  precipirat ion  module  is  called  and  it  determines 
the  intensity  (qualitatively)  of  the  precipitation. 

The  cloud  top  temperature  module,  in  which  the 
temperature  and  the  height  of  the  clouds  are  calculated,  is 
called  for  all  cloud  cases.  The  FNOC  surface  temperature 
value  is  used  in  the  cloud  top  temperature  calculation, 
while  t he  FNOC  upper  air  data  are  used  in  the  height 


calculation 


TABLE  XIII 

SPADS  Cloud  Hodel  Precipitation  Intensity  Categories 

Equation  sual  ■  IR  radiance  +  average  VIS  brightness 

No  rain  ($)  suml  <  184 

Light  rain  (1)  184  <  suml  £  195 

Moderate  rain  (2)  195  <  suml  <_  224 

Heavy  rain  (3)  224  <  suml 

The  final  step  is  processing  the  results.  cloud  types 
(Table  XII),  precipitation  intensities  (Table  XIII),  cloud 
top  temperatures  (K)  ,  and  cloud  top  heights  (mb)  are  stored 
in  a  permanent  file  for  future  display.  The  average  visual 
brightness,  standard  deviation,  and  the  amount  of  cloud  that 
correspond  to  each  IR  pixel  are  printed  for  reference  and 
validation  purposes. 
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A.  INTRODUCTION 

A  test  was  conducted  to  dateraine  if  the  SPADS  Cloud 
Model  concepts  were  valid.  The  test  was  designed:  (1)  to 
capture  a  small  area  cf  coincident  IR  and  vis  satellite  data 
in  the  mid-latitudes  and  to  obtain  FNQC  fields  for  the  same 
date-time  group  and  location  as  the  satellite  data,  (2)  to 
run  the  SPADS  Cloud  Model  using  these  data,  and  (3)  to  ana¬ 
lyze  the  results  of  the  model  to  determine  the  validity  of 
the  model* s  analyses.  The  following  sections  describe  the 
input  to  this  test,  test  procedures,  results,  and  subsequent 

9 

iterations  and  adjustments  made  to  improve  the  model's 
performance. 

B.  DATA 

GOES  IR  and  VIS  data  were  captured  and  archived  on 
November  9,  1982  at  1315  local  time,  for  a  site  centered  at 
location  latitude  44° N,  longitude  14  1°W.  The  corresponding 
FNOC  data  fields  were  also  saved  for  the  same  date  from  1200 
GMT  data.  Fig.  34  is  the  GOES  IR  image  from  which  an  area  of 
approximately  460  X  460  n  mi.  was  selected  for  analysis. 


This  afforded  4096  IR  pixels  for  evaluation.  Fig.  35  is  the 


IR  satellite  picture  of  the  selected  area  with  the  center 
point  matching  that  in  ?ig.  34.  The  IR  data  resolution  was 
approximately  7  n  ai.  at  this  latitude  and  4  n  ai.  at  the 
subsatellite  point.  The  visual  data  resolution  used  was  0.9 
n  ai.  (0.5  n  ni.  at  subsatellite  point)  .  The  visual  area 
(see  Fig.  36)  is  the  sane  area  as  the  infrared  giving  an  8X8 
visual  grid  of  pixels  (64  points)  per  IR  pixel. 

The  FNOC  vertical  temperature  profile  was  obtained  to 
match  the  center  point  of  the  satellite  pictures  and  the 
corresponding  heights  (mb)  are  listed  in  Table  XIV.  Only  one 
set  cf  temperatures  was  used  because  the  resolution  of  the 
FNOC  fields  (63  X  6  3  field  with  381  ki  mesh  length)  is 
approximately  the  same  as  the  area  of  the  satellite 
pictures  used  in  this  test. 

C.  TEST  PROCEDURES 

The  following  list  is  the  set  of  procedures  for  this 
test  of  the  SPADS  Cloud  Model.  The  steps  taken  were  to: 

1.  Capture  coincident  VIS  and  IR  satellite  images. 

2.  Receive  surface  and  upper  air  r.emoaratures  from 
FNOC. 

3.  Manually  analyze  images  for  cloud  boundaries  and 
types. 


4.  Select  the  temperature  profile  closest  to  satellite 


5.  Ran  SPADS  cloud  Model  computer  program. 

6.  Output  results  in  a  fora  at  f.ot  analysis. 

7.  Coapare  manual  tc  computer  output. 

8.  Adjust  computer  model  as  appropriate. 

9.  Repeat  steps  5,  6,  and  7. 

The  manual  nephanalysis  was  completed  by  Dr.  C.H.  Wash, 
Assistant  Professor  in  Meteorology  at  the  Naval  Postgraduate 
School.  He  used  both  the  IR  and  Visual  images  (Fig.  34,  35, 
and  36)  to  discern  the  cloud  boundaries  and  types  seen  in 
Fig.  37  (see  Table  XV  for  symbol  definition) 

During  the  first  computer  run,  a  cloud  threshold  of  5 
visual  digital  counts  was  used.  This  value  corresponded  to 
the  lowest  visual  threshold  that  could  have  clouds  (see  Fig. 
33)  .  After  analysis  of  the  first  run,  however,  a  threshold 
of  20  counts  was  selected.  This  adjustment  accounts  for  our 
uncertainty  in  the  GCES  VISSR  visible  channel  calibration, 
as  compared  to  that  for  the  TIROS  AVHRR  visible  channel.  In 
normal  use,  VISSR  is  a  relatively  qualitative  radiometer, 
and  the  engineering  data  needed  to  convert  values  to  high- 
quality  calibrated  radiances  were  not  readily  available. 
Two  sets  of  standard  deviations  were  used  and  are  presented 
in  Table  XVI. 
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TABLE  ZZY 

F»oc  T..p.r«t«jet||o|ilj(|tjr  9<  1,82  at 


Standard  level  (mb) 
100 
150 
200 
250 
300 
400 
500 
700 
850 
1000 

1010  (SFC) 


Temperature  (°K) 

212. 

214. 

216.8 

221.8 

229.7 

244.6 

254.3 

267.7 
275.5 

283.4 

289.7 
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TAB LB  XV 

Definitions  of  Syebols  Deed  in  Fig*  37 


Nimbostratus  -  Multiple  layers 

NS -ML 

A1 tostratus 

AS 

Altocunul  us 

AC 

Cirrus 

Ci 

Cunul  us 

Cu 

Stratus 

St 

Stratocunulus 

SC 

Broken 

© 

TABLE  XVI 

Standard  Deviation  (SXGi)  Values.  Used  for  Run  1  and  2  (see 
Table  XII  for  how  they  are  used) 

SIG  (1)  Teat  1 _ Test  2 

4  2.5 

20  20 
3  1.5 


SIG  1 
SIG  2 
SIG  3 


D.  RESULTS 

The  following  subsections  discuss  the  results  of  the  t 
computer  runs  as  compared  to  the  manual  analysis.  Th9 


results  are  discussed  by  the  type  of  product  and  run  number; 


where  the  products  are  cloud  type,  precipitation,  cloud 
aaount,  cloud  top  teeperature  and  height,  and  texture 
(standard  delation)  . 

1,  Run  2 

a.  Cloud  Typing 

The  overall  general  patterns  of  the  cloud  types 
(see  Fig.  38)  cospared  well  to  the  eanual  analysis  (see  Fig. 
37).  The  observed  cloud  structure  is  most  easily  discussed 
by  considering  the  inages  as  four  quadrants;  northeast  (NE)  , 
northwest  (NW)  ,  southwest  (SJf)  ,  and  southeast  (SE)  with 
north  being  the  Top  of  the  iaag9 . 

The  NE  quadrant  has  aid-level  clouds  such  as 
altostratus  (see  Fig.  37).  The  SPADS  cloud  aodel  did  well 
here  indicating  predoainatly  aid-level  clouds  (see  Fig.  38)  . 
The  aodel  also  aatched  the  aanual  analysis  in  the  3W  quad¬ 
rant.  This  area  has  high  cloud  tops  of  niabostratus  and 
embedded  cuauloniabus . 

The  SH  quadrant  has  a  mixture  of  aid-level  and 
aulti- layered  thick  clouds  and  again  the  cloud  aodel  is  in 
agreeaent  with  the  aanual  analysis  (see  Fig.  37  and  38).  The 
objective  analysis  shows  sore  detail  than  the  aanual 
analysis  because  the  coaputer  program  is  working  pixel  by 


pixel  while  the  manual  analysis  can  only  resolve  groups  of 
pixel  defining  cloud  areas. 

In  the  SE  quadrant,  the  SPADS  cloud  aodel  could 
not  resolve  the  difference  between  low-level  (stratus)  and 
aid- le«el  clouds.  The  boundary  between  these  is  located  on 
the  upper  border  of  the  SB'  quadrant.  Another  discrimination 
problem  is  between  stratus  and  cuaulus  in  the  lower  right 
quadrant.  The  amount  of  cuaulus  humulus  that  seeas  apparent 
in  the  visual  satellite  picture  (see  Fig.  36)  and  the  annual 
analysis  (see  Fig.  37)  did  not  appear  in  the  coaputer 
analysis  (see  Fig.  38).  This  problem  is  related  to  the  tex¬ 
ture  threshold  tests,  as  borh  stratus  and  cumulus  humulus 
have  the  same  cloud  type  thresholds, 
b.  Precipitation 

The  patterns  of  rain  intensity  (see  Fig.  39)  in 
the  NW  and  SW  quadrants  are  based  on  the  premise  that  the 
thicker  and  colder,  therefore  brighter,  clouds  (compare  Fig. 
39  to  Pig.  35  and  36)  have  a  greater  probability  of  rain  and 
a  higher  intensity  of  rain.  Along  with  this,  is  the  knowl¬ 
edge  that  the  approximate  threshold  for  rain  clouds  is 
colder  than  -22  °C  (Liljas,  1981a;  Barrett  and  Martin, 
1981).  3oth  coaputer  runs  gave  the  saae  agreement  with  the 
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Figure  38.  Cloud  Types  as  Analyzed  by  ths  SPADS  Cloud  Scdel 
in  Sun  i  (see  Table  XII  for  numeric 
definitions)  . 


manual  analysis  bee  a use  the  areas  of  rain  clouds  did  not 
change.  it  this  time,  no  surface  data  are  available  for 
ground-truth  verification  of  the  rain  patterns.  “ 

c.  Cloud  Aaount 

The  cloud  asount  for  the  first  run  was  100%  for 

I 

the  entire  analysis  area.  The  visual,  IB,  and  annual  pic¬ 
tures  all  shewed  clear  and  scattered  cloud  areas  in  the  SE 
quadrant.  In  the  annual  analysis  (see  Pig.  37),  this  quad¬ 
rant  has  these  sections  marked  Cu,  St/Sc,  and  o  is  (see 
Table  XV  for  symbol  definitions)  .  Further  research  into  the 
cloud  threshold  for  GOES  satellite  iaagery,  combined  with 
analysis  of  the  high  resolution  visual  data  used  in  this 
test,  suggested  a  better  cloud  threshold  value  would  be  20 
visual  counts  (Huench  and  Keegan,  1979,  and  see  discussion 
in  Section  C) .  The  results  from  using  this  threshold  are 
discussed  in  the  Hun  2  section. 

d.  Cloud  Top  Temperatures  and  Heights 

The  cloud  top  temperatures  (see  Pig.  40)  and 
heights  (see  Pig.  41)  look  reasonable  for  the  cloud  types 
identified.  The  HE  quadrant  had  temperatures  ranging  from 
260  to  250  X  and  a  height  of  a pproixmately  500  mb,  these 
matching  the  type  of  mid-level  cloud.  The  NW  guadrant 
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clouds  are  mostly  nimfcostratus  with  high-level,  cold  cloud 
tops.  The  height  is  appro ximately  300  ab  and  the 
temperature  is  230  K.  Again,  agreeaent  is  good.  The  aired 
clouds  in  the  SW  quadrant  also  have  aired  height  (aid-level 
and  high-level)  and  teaperatures  which  coapare  well  (see 
Fig.  40  and  41) 

The  SE  quadrant's  heights,  ranging  from  850  to 
700  mb,  define  the  area  reasonably  well  (see  Fig.  41).  The 
division  for  the  low  and  aid-level  clouds  (see  Pig.  37) 
appears  as  the  700  ab  contour  line  (see  Fig.  40).  The  cloud 
teaperature  field  (see  Fig.  40)  defines  the  division  of  low 
and  mid- level  clouds  less  clearly,  but  partly  cloudy  areas 
(275  K  contour)  are  indicated.  The  cirrus  and  altostratus 
clouds  are  too  cold  and  high,  as  compared  to  the  FNOC 
temperature  profile.  This  error  appears  to  be  in  the  appli¬ 
cation  of  the  bispectral  radiance  equation  (Eq.  4,  Table  IV) 
from  Reynolds  and  Vender  Haar  (1979)  in  partly  cloudy 
regions.  The  coefficients  in  that  equation  depend  strongly 
on  the  radiance  calibration  of  the  radiometer  used.  The 
assumed  calibration  coefficeints  used  in  this  study  are 
clearly  inadequate  and  must  be  improved  through  empirical 
validation  of  GOES  VISSR  calibration  using  ground  truth 


Figure  41.  Cloud  Top  Heights  (mb)  from  Run  1  of  the  SPADS 
Cloud  Model. 


e.  cloud  Texture  (Standard  Deviation) 


The  cloud  texture  (standard  'deviation) 
calculation  results  are  presented  in  Fig.  42.  Overall,  they 
show  an  expected  pattern;  nixed  and  broken  cloud  areas  (in 
the  SI  and  SE  quadrants)  have  the  highest  variation  and  uni** 
form  stratus  type  cloud  cover  areas  (in  the  HI  and  HE  quad¬ 
rants)  have  the  smallest  variation  (coepare  Fig.  36  and  42). 
The  highest  deviations  (approximately  4.0)  were  located  in 
the  SE  guadrant  in  the  areas  discussed  previously  under 
cloud  amount,  after  contouring,  there  appeared  to  be  a  defi¬ 
nite  boundary  between  the  highly  variable  region  (SE  quad¬ 
rant)  and  the  nearly  uniform  areas  (HE  quadrant,  see  Fig. 
42) .  The  SE  quadrant  had  been  manually  typed  as  cumulus, 
with  broken  areas,  and  stratus/stratocumulus  clouds  (see 
Fiq.  37)  .  After  evaluation,  the  standard  deviation  threshold 
value  for  discerning  cumulus  humulus  was  adjusted  along  with 
the  value  for  altostratus.  There  did  not  appear  to  be  obvi¬ 
ous  patterns  of  different  size  cumulus  congestus  and  there¬ 
fore  this  standard  deviation  threshold  was  not  adjusted. 
The  standard  deviation  definitely  shows  promise  in 


discriminating  between  cumuloform  and  stratus  cloud  types 


Cloud  Texture  (Standard  Deviation  of  the  Digital 
Counts)  from  Run  1  of  the  SPADS  Cloud  Model. 


The  results  from  the  second  run  are  presented  below 
with  the  revised  cloud  and  standard  deviation  thresholds. 
Only  regions  where  changes  from  Run  1  occurred  are  discussed 
in  detail. 

a.  cloud  Typing 

The  cloud  typing  results  were  changed  soaewhat 
because  of  the  new  cloud  threshold  and  standard  deviation 
thresholds  (see  Fig.  43)  especially  in  the  lower  right  quad¬ 
rant.  More  cumulus  huaulus  and  altostatus  types  were  deli- 
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neated  and  the  patterns  appear  closer  to  those  in  the  manual 
analysis. 

b.  Cloud  Amount 

The  cloud  amount  patterns  (see  Fig.  44)  are  more 
in  line  with  the  manual  analysis  and  what  can  clearly  be 
seen  on  the  satellite  images.  It  appears  thatthe  cloud 
threshold  of  20  is  better,  but  this  may  require  additional 
study. 

c.  Cloud  Top  Temperature  and  Height 

The  problem  in  obtaining  the  correct  cloud  top 
temperature  and  height  in  a  partly  cloudy  situation  using  a 
bispectral  technique  still  remained.  This  problem  stems 
from  uncertainty  in  calibration  of  the  GOES  7ISSR.  The  VISS3 
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Figure  4  3. 


Cloud  Types  as 
Model,  hn  2  { 
definitions)  . 


Figure  44.  Cloud  Amount  (in  lOOths)  from  Run  2  of  the  SPADS 
Cloud  Model  (the  number  in  the  enclosed  areas 
represents  the  lowest  cloud  amount  for  the  area) . 


infrared  channel  is  not  actively  calibrated  and  extensive 
ground  truth  is  required  to  accurately  interpret  its  data  in 
teres  of  teeperatures  (Haul#  198  1).  In  the  present  progran# 
teaperature  calibration  would  an  ter  into  the  coefficients  of 
the  cloud  height  Eg.  4  (Table  17)  and  directly  into  the 
cloud  top  teeperatures.  Tiae  was  not  available  in  this  study 


to  make  this  adjustaent 
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A.  DISCOSSICH 

This  chapter  saaaari zes  the  thesis  accomplishments,  the 
problees  encountered,  and  suggestions  for  future  research. 

The  first  accoaplishaent  was  to  complete  an  overview 
of  cloud  and  precipitation  studies  with  an  eephasis  given  to 
those  using  an  interactive  computer  system.  Next,  useful 
techniques  and  ideas  were  selected  froa  these  studies  and 
used  to  design  the  SPADS  cloud  and  precipitation  aodel.  This 
SPAOS  Cloud  Hodel  was  then  coded  and  established  on  the 
SPADS  system.  Finally,  two  test  runs  on  a  test  case  were 
completed  with  good  success. 

Within  an  hour,  detailed  information  was  produced  by 
the  SPADS  Cloud  Model  through  concurrent  analysis  of  high 
resolution  visual  and  infrared  digital  satellite  data.  Over 
a  maritime  region,  the  SPADS  Cloud  Model  estimated  various 
cloud  types,  precipitation  intensity,  cloud  amounts,  and 
cloud  top  temperatures  and  heights. 


2.  gSgfelftSg 


Several  probleas  were  encountered.  The  cloud  mount 
threshold  needed  adjustaent  and  two  tests  of  different 
thresholds  were  acccaplished.  The  threshold  of  20  was 
satisfactory  but  further  verification  is  suggested. 

Another  area  cf  concern  is  the  threshold  values  for 
the  texture  tests.  Two  thresholds  were  adjusted  ir.  the  tests 
but  all  three  need  further  refineaent.  Ihe  separation  test 
(the  least  square  fit  equations)  for-  cirrus  and  altostratus 
could  not  be  substantiated  in  this  data  set  and  will  also 
need  verification  in  future  runs. 

An  additional  problem  is  in  the  cloud  top 
temperature  and  height  calculations.  As  discussed  before, 
the  temperatures  are  too  cold  and  the  heights  too  high  in 
certain  partly  cloudy  areas.  The  bispectral  equation  to 
obtain  these  values  needs  further  study. 

All  three  of  these  probleas  are  related,  at  least  in 
part,  to  the  problem  of  inter-calibration  of  AVHRR  and  GOES 
VISSR  visible  and  infrared  data. 
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The  following  are  reco eiendatioas  and  suggestions 


for  additional  research.  The  SPA  OS  Cloud  Model  should  be: 


a.  converted  to  use  temperatures  for  infrared  thresholds 
and  albedoes  for  visual  thresholds  to  allow  easier  con* 
version  between  different  satellites. (polar  and  geosta¬ 
tionary)  and  seasons  (this  will  require  developaent  and 
addition  of  sensor  calibration  tables  and  nodules  to 
the  SPADS  Cloud  Model). 


b.  converted  tc  a  variable  resolution  so  that  larger  areas 
could  be  analyzed  in  the  sane  tine  fraae  but  also  keep 
the  high  resolution  capability. 


c.  tested  and  verified  m  various  regions  by  using  surface 
observations,  rawinsondes,  aircraft  platform 
observations,  ship  observations,  radar,  and, if  availa¬ 
ble,  the  0.  S.  Air  Force  3DNEPH  cloud  analysis. 


Digital  satellite  data  are  currently  underutilized. 


except  perhaps  for  sea  surface  temperature  and  temperature 


sounding  applications.  Most  visual  and  infrared  satellite 


data  are  treated  as  images  and  are  still  analyzed  manually. 


Now,  with  the  development  of  SPADS  and  similiar  systems,  new 


satellite  products  can  be  produced  in  regional  forecast 


sites  based  on  digital  satellite  data.  This  thesis  illus 


trates  the  detailed  analysis  of  the  digital  satellite  infor¬ 


mation  that  can  be  received  by  the  user  in  a  short  (1-2 


hours)  time  fraae.  Specifically,  this  work  was  directed  to 


using  the  SPADS  computer  system  toward  process  satellite 


visual  and  infrared  digital  data  and  to  produce  cloud  and 


precipitation  information.  As  SPADS  is  placed  in  the  field 
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forecast  sites,  this  information  will,  for  the  first  tine, 
be  eade  available  in  an  operational  setting  instead  of  being 
available  only  to  aeteorologists  at  the  aajor  satellite  cen¬ 
ters.  In  the  future,  this  and  other  applications  will  sub¬ 
stantially  aid  the  operational  forecaster. 
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APPENDIX  A 

SPADS  CLOUD  1SODEL  COMPUTER  PROGRAM 


C  CLOD  —  IS  DRIVER/?1AIN  PROGRAM 

C  THIS  PROGRAM  ANALYZES  VIS  AND  IP.  TOGETHER  TO  DERIVE  TYPES  OF 
C  CLOUDS, CLOUD  HEIGHTS , CLOUD  AMOUNT, CLOUD  TOP  TEMPERATURE, 

C  AND  RELATIVE  INTENSITY  OF  RAIN. 

C  PROGRAMMED  3Y  LT.  C.  NELSON 

C  INITIALIZATION 

C 

COMMON/ A/ML(128)  .CIRRUS,  1C, KI.ISEL 
COMMON/ 3/ CTYPE ( 64 ) ,PRCIPI(64) ,MSAVG(64) ,M 
COMMON/ C/SIGMA( 64) ,SIGI , SIG2 ,SIG3 

COMMON/ D/ HEIGHT (64) ,CTEMP(64) ,ACDAVG(64) ,E1 , E2 , ALAT , ALON , J 1 , 
*TLV(11),E3 

C0MM0N/F/ZP,XCEN,YCEN,FLIP,A1 ,3ASL0N,LVL( 11) 

COMMON/ 0/ STOP (11) 

REAL  ML,MS2(8,S) ,MSAVG, SIGMA, SI ,52,5, ACLD,ACnAVG,AMCLD, SUM, SUMSIG, 
♦HEIGHT , N , CTF.MP , S IG l ,  SIG2 , SIG3 

INTEGER  CTYPE, CIRRUS, PRCIPI ,M,0,P ,L,Dr,LDl ,PCLD2 ,0,  TSF,L,MS(  256 ,") , 
*VISFN(7) ,IRFN(7) ,MSl (3) 

r 

11*1 
ISEL-0 
''1=0 
'1=0.0 
El =0.9 
E2=0. 55 
E3»0.0 
sum-o.o 
SUMSIG-0.0 
AMCLD-0,0 
CIRRUS-0 
ACLD-0 . 0 
?=64 

SIOl-2.5 

SIG2-20.0 

SIG3-1.5 

PCLDL-20.0 

0=512 


« 


non 


MENU 


-J*  _*  — *  -A,-.  JL.  -  - 


P 


TYPE  "<15>" 

TYPE  "THIS  PROGRAM  IDENTIFIES  CLOUD  TYPES" 

TYPE  "ADDITIONAL  CLOUD  DATA  CAN  ALSO  BE  OBTAINED" 
TYPE  "THESE  ARE  -  " 

TYPE  "0  NO  ADDITIONAL  INFORMATION" 

TYPE  "I  PRECIPITATION  INTENSITY" 

TYPE  "2  CLOUD  TOP  TEMPERATURE" 

TYPE  "3  CLOUD  TOP  HEIGHT" 

TYPE  "4  ALL  OF  THE  ABOVE" 

TYPE  "3  EXIT  PROGRAM" 

ACCEPT  "ENTER  REQUEST-  ",ISF,L 
IF(ISEL.EQ.8)G0  TO  80 
IF( I3EL.EQ.0.0R. ISEL.EQ. 1 )G0  TO  131 
CALL  0PEN(4 ,"TLVN.O" , 1 ,IER) 

IF( IER.EQ. 1 )G0  TO  170 
TYPE  "TLV  FILE  OPEN  ERROR,  RC-  ",IER 
170  CALL  RD8LK(4 , 1 ,MS , 1 , IER) 

IF( IER.EQ. l)GO  TO  ISO 
TYPE  "TLV  FILE  READ  ERROR,  RC-  ",IER 
GO  TO  70 
ISO  DO  l  1-1,11 

TLV(I)»FLOAT(MS(I))/100.0 

1  CONTINUE 

VRITE( 1 2 , 826) ( TLV( L) ,L- 1 , l 1 ) 

826  FORMAT ( 2X.11F6.2) 

2  FORMAT( S13) 

131  TYPE  "<15>" 

ACCEPT  "ENTER  IR  FILENAME" 

RF.AD(  1 1 ,2)  IRFN(  1 ) 

CALL  OPEN(2,IRFN,l,IER) 

IF( IER.EQ. l)GO  TO  IDO 
TYPE  "IR  PILE  -  OPEN  ERROR,  RC-  " , IER 
ACCEPT  "TRY  AGAIN?", IX 
IF( IX. EQ. 1)00  TO  IS  1 
IF( IX. SO. 0)00  TO  70 
IDO  ACCEPT  "ENTER  VIS  FILENAME" 

READ( 1 l , 2) vI8rN( 1 ) 


123 


CALL  OPEN( 1 ,VISFN, 1 ,IER) 

IF(IER.E0.1)GO  TO  200 

TYPE  "VIS  FILE  -  OPEN  ERROR,  RC-  ",IER 

ACCEPT  "TRY  ACAIN?",IX 

IF( IX. EO. 1 )G0  TO  190 

GO  TO  69 

OPEN  FILES  FOR  STORAGE  OF  CALCULATED  VALUES 

200  CALL  0PEN(5 ,"CALC1D.0" , 3 , IER) 

IF( IER.EO. 1)G0  TO  201 

TT°P.  "0AT\  1  FILE  OPT'  ERROR,  RC-  ",IER 
GO  TO  6R 

201  CALL  0?EN(7,"CALC2D.Q",3,IER) 

IF( IER.EQ. 1)G0  TO  202 

TYPE  "DATA  2  FILE  OPEN  ERROR,  RC-  ",I?.R 
GO  TO  69 

DO  LOOP  TO  CALL  IN  SERIES  OF  DATA 
".LOCKS  RUN  FROM  0  TO  511 

LOOP  STARTS  AT  192,192  FOR  IR  TO  MATCH  VIS  (COLLOCATED  CENTER 
POINT).  IR  IS  2  MILE  THEREFORE  NEED  TO  USE  EVERY  OTHER  VALUE 

202  13-0 
12-0 

DO  32  11-1 ,120,2 
I3-(  Il-D+191 
I2-1+I2 
DO  6  I-l,P 
HEIGHT( I)«0.0 
PRCIPI( I)-0.0 
ACDAVG(I)-0.0 
SIGMA( I)-0.0 
CTYPE(I)-0 
CTE’D’C  I)-0.0 
6  CONTINUE 

DRIVER  MODULE  WHICH  WILL  PREPARE  VISUAL  DATA  TO  RE  USED  IN  OTHER 
MODULES.  THE  AVERAGE  VISUAL  BRIGHTNESS,  THE  STAWARD  DEVIATION  OF 
VISUAL  BRIGHTNESS,  .AND  THE  CLOUD  AMONT  IN  R”S  GRID  SETS  ''ILL  3E 


DETERMINE. 

IR  DIGITAL  VALUE  FILE 


CALL  RDBLK(2,I3,MS(l,l),l,IER) 

IF( IER.EQ. 1)G0  TO  203 

TYPE  "IR  FILE  -  READ  ERROR,  RC-  ",IF.R 

GO  TO  69 

203  CALL  IJPK(MS(96 ,1)  ,MS(  1,2), 128) 

DO  100  1-1,123 
M.L(I)-FLOAT(US(I,2)) 

100  CONTINUE 
XI -0 

Jl-(I2-l)*3 

CALL  RD3L’<(  1  ,J1  ,MS,3 ,IER) 

I?( IER.EO. l)GO  TO  204 

TYPE  "VISUAL  FILE  READ  ERROR,  RC-  ",IER 
GO  TO  69 

204  CONTINUE 

DO  30  '<2-1,128,2 

Kl-Kl+1 

X-K2+1 

rilS  LOOP  PROCESSES  VISUAL  DATA  3Y  SX3  OR  ID 


DO  3  1-1,3 
US  1 ( I ) -0 . 0 
DO  7  J2-1 ,3 
US2(  I, J2J-0.0 
7  CONTINUE 
3  CONTINUE 
Ll«(Kl-l)*4+l 
DO  11  0-1,8 

CALL  UPX(MS(L1  ,0)  ,  IS  1 ,3) 


DO  10  L-1,3 

U.S2(L,0)-FL0AT(  ISl(L)) 

10  CONTINUE 

11  CONTINUE 

Nl-8 

DO  13  L-1,3 
do  i:  o-l  *;i 


N-Nfl.O 

SUM-SUM+MS2(L,Q) 

SUMSIG-SUMSIG+(MS2(L,Q))**2.0 

IF(MS2(L,O).GE.?CLDl)ACLD-1.0 

IF(MS2(L,Q).LT.PCLJ)l)ACLD-0.0 

AMCLD-ACLtH-AMCLD 

12  CONTINUE 

13  CONTINUE 

SIOMA( K1 )»SQRT( (N*SUMSIG-SUM**2 .0) /( M*( N-l .0) ) ) 

MSAVG(K1)«SUM/N 

CALL  CLOUO 

ACDAVG(K1  )-AHCLD/!I 

IK\CTYPE(  K1 )  .  SC<.3.0R.CTYPE(K1 )  .EO.  9)  £3-1 .0 
I?(ISEL.EQ.O.OR.ISEL.RQ.1)GO  TO  13 
CALL  TRMMT 

IF(ISEL.EQ.3)GO  TO  69 

13  cmus-o 

N-0.0 
ACLD-0.0 
AMCLO-O.O 
SUM-0. 0 
SUMSIG-0.0 
E3-0.0 
30  CONTINUE 


STORE  ALL  VALUES 

VALUES  ARE  STORED  IN  CALCl  PILE  —  CTYP* ,°EClPt ,CTEMP , 
HEIG3TC 1-256  VD) 

VALUES  IN  CALC2  ~  HSAVC, SIGMA, ACDA7G,0( 1-256  WD) 


DO  23  1-1,64 
MS( I, 1)-CTYPE( I) 

J-I+64 

MS( J, 1)-PRCIPI( I) 

T5-J+64 

MS(<5,l)-I?IX(CTEMP(I)+.5) 

K4-K5+64 

MS(K4,l)«IFIX(-lEIOWT(I)+.5) 

MS( I , 2 ) -IF IX( MS A VC ( I )+ . 5 ) 
MS(J,2)-IFIX(100.0*SICMA(I)+.5) 


m 


MS(K5,2)«IFIX(100.0*ACDAVG(I)+.5) 

MS(K4,2)-0 

23  CONTINUE 
K7-I2-1 

CALL  WRBLK(5,K7,MS(ltI),l,IER) 
IF(IER.EQ.l)GO  TO  22 
TYPE  "WRITE  ERROR  5  RE-  "f IER 
GO  TO  69 

22  CALL  WRBLK(7,K7,MS<1,2),1,IER) 

IF( IER.EO. 1 )GO  TO  26 

TYPE  "WRITE  ERROR  7  RC-  ",IER 

GO  TO  60 

26  WRITE (  12, 27)('-tEICHT(I),  1-1,64) 
WRITE ( 1 2 , 24 ) ( S IGWA( I ) , I- 1 , 64 ) 
WRITE ( 1 2 , 27 ) ( CTEMP( I ) , I- 1 , 64 ) 
WRITE( 1 2 , 24 ) ( MSAVG( I ) , I- 1 , 64 ) 
WRITE ( 12 ,24) ( ACDAVG( I) ,1-1,64) 
WRITE( 12 ,25) (CTY°E( I) ,1-1 ,64) 
WRITE( 12 ,25) ( PRCI?I( I) ,1-1,64) 
WRITE(  12, 27 )(ML( I) ,1*2,123,2) 

24  FORMAT(/,4(/,I6(IX,F7.4))) 

25  FORMAT( / ,6412) 

27  FORMAT( /,4(/,16(lX,F7.2))) 

ACCEPT  "CONTINUE?" , 1ST? 

I?( 1ST?. EO, I )GO  TO  69 

32' CONTINUE 

69  CALL  CLOSE(l.IER) 

70  CALL  CLOSE( 2 , IER) 

CALL  CLOSE (4, IER) 

CALL  CLOSE(5 ,IER) 

CALL  CLOSE(7 ,IER) 

30  STOP 


a 


k.  -V. 


_ _ _ _ _ . 


SUBROUTINE  CLOUD 

CLOUD  TYPING  MODULE,  ALSO  THE  PRECIPITATION  MODULE  IS  CALLED 
FROM  HERE  WHEN  PRECIPITATION  CLOUDS  ARE  IDENTIFIED. 

COMMON/ A/ML( 128), CIRRUS, K.K1.ISEL 
COMMON/B/CTYPE(64) ,PRCIPI(64) ,MSAVG(64) ,M 
REAL  ML ,MSAVG , HEIGHT , CTEMP , SIG , SIG 1 , SIG2 , SIG3 
INTEGER  I , CTYPE ,PRCIPI , CIRRUS , ISEL 


THESE  ARE  LILJAS  CONVERTED  GOES  THRESHOLDS 


H2-  *4.0 
H3-100.0 
H4-U8.0 
H5-159.0 
H6-255.0 

THE  FOLLOWING  IF  STATEMENTS  CHECK  BOTH  IR  (ML)  AND  VIS  (MSAVG) 
VALUES  AGAINST  THRESHOLDS  TO  DETERMINE  CLOUD  TYPE.  THE 
STANDARD  DEVIATION  (SIG)  OF  THE  3X8  VISUAL  GRID  IS  USED. 

TYPE  1  AND  2  -  CIRRUS/ ALTOSTRATU S 

IF(ML(’<)  .GT.H4.AND.ML(K)  .LE.H6.AND.MSAVG(K1 )  .GT. 5.0. AND. 
*MSAVG(K1).LE.23.0)GO  TO  700 

TYPE  3 - FOG/ STRATUS /CUmUS  HUMIJLUS 

IF(ML(K). CT.H2. AND. ML(K).LE,H4.. AND. MSAVG(K1).GT. 13.0. AMD. 
*MSAVG(K1).LE.20.0)GO  TO  900 

TYPE  4  AND  5 - THICK  FOG /STRATOCIDIULUS/ CUMULUS  UUTTJLUS 

IF(ML(K) .GT.H3.AND.ML(K) .LE.H4 . AND.MSAVG(K1 ) ,G7. 20.0. AND. 
*MSAVG(K1).LE.36.0)GO  TO  600 
GO  TO  601 

600  M«1 

GO  TO  800 

TYPE  6  AND  7 - SMALL/LARGE  CUMULUS  CONGESTTJS 

601  IF(ML(K)  .GT.H4.AND..'IL(K)  .LE.H5.A?TD.MSAVG(K1 )  .GT.  23.0.  AND. 
*MSAVG(K1).LE.39.0)G0  TO  605 

GO  TO  606 
605  M-2 

GO  TO  800 


a  o 


C  TYPE  8  -  NIMBOSTRATUS  (MULTI-LAYERED) 

606  IF(ML(K) .GT.H5.AND.ML(E) .LE.H6.AND.MSAVG(K1) .GT. 23.0. AND. 
*MSAVG(K1).LE.31.0)GO  TO  910 

C  TYPE  9  -  CUMULONIMBUS 

IF(ML(K).GT.H5.AND.ML(K).LE.H6.AND.MSAVG(K1).GT.31.0) 

*G0  TO  915 
CTYPE(Kl)-0 
GO  TO  38 

C  CALL  STDDEV  MODULE  TO  FURTHER  DISCRIMINATE  CLOUD  TYPES 
700  M-3 

800  CALL  STDDEV 
GO  TO  38 

C  LABEL  TYPES 
900  M-4 

GO  TO  800 
38  PRCIPI(Kl)-0 
GO  TO  40 
910  CTY?E(Kl)-8 
GO  TO  920 
915  CTYPE(Kl)-9 

BOTH  8  AND  9  TYPES  ARE  CONSIDERED  RAIN  CLOUDS,  SO  CALL 
PRECIP  MODULE. 

920  IF( ISEL.EO.0.OR.ISEL.EQ.2.OR.ISEL.EQ. 3)G0  TO  40 
CALL  PRECIP 
40  RETURN 
END 


SUBROUTINE  STDDEV 


STANDARD  DEVIATION  TEST  MODULE  HELPS  DISCRIMINATE  FURTHER  THE 
TYPE  OF  CLOUDS.  ALSO  IT  SETS  CIRRUS  FLAG  FOR  USE  IN  THE  TEHHT 
MODULE. 

COMMON/ A/ML( 128) .CIRRUS ,K,K1 ,ISEL 
COMMON/B/CTYPE(64) ,PRCIPI(64) ,MSAVG(64) ,M 
COMMON/ C/SIGMA(64) ,SIG1 ,SIG2 ,SIG3 
REAL  ML , MSAVG , HEIGHT , CTEMP , SIG , Y , SIG 1 , SIG2 , SIG3 
INTEGER  CTYPE.PRCIPI, CIRRUS, I, M 


IF(M.EQ.4)GO  TO  SIS 
IF(M.EQ. l)GO  TO  810 
IF(M.EO. 2)GO  TO  813 
CIRRUS  VS  ALTOSTRATUS,  FIRST  TEST 

THIS  EQUATION  IS  FOR  LIUAS  CONVERTED  GOES  THRESHOLDS 
Y- 1 58 . 094-6 . 807  6*MSAVG( K1 )+4 . 1 59*MSA VG( El )**2 . 0 
IF(Y.LE.ML(K) )G0  TO  814 
SECOND  TEST  FOR  CI/AS 

THIS  EQUATION  IS  FOR  LILJAS  CONVERTED  GOES  THRESHOLDS 

Y- 1 1 5 . 24+0 . 5246*MSAVG( K1 )+0 . 0395*MSAVG( XI ) **  2 . 0 

IF(Y.CE.ML(K))CO  TO  816 

T1IRD  TEST  FOR  CI/AS 

IF(SIGMA(K1).GT.SIG3)G0  TO  314 

CO  TO  816 

810  I?(SIGMA(K1).LT.SIG1)G0  TO  811 
TYPE  5  —  CUMULUS  HUMULUS 
CTYPE(Kl)-5 
GO  TO  820 

TY»E  4  —  THICK  FOG/STRATOCUMULUS 
311  CTYPE(K1)«4 
GO  TO  820 

313  IF(SIGMA(X1).LT.SIG2)G0  TO  817 
CTYPE(Kl)-6 

TYPE  6  —  SMALL  CUMULUS  CONGESTiJS 

GO  TO  820 

TYPE  l  —  CIRRUS 

314  CTYPE(K1)-1 
CIRRUS- 1 


38 


GO  TO  820 

TYPE  2  —  ALTOSTRATUS 

816  CTYPE(Kl)-2 
GO  TO  820 

TYPE  7  —  CUMULUS  CONGE STUS  - 

817  CTYPE0Cl)-7 
GO  TO  820 

818  IF(SIGMA(K1).LT.SIG1)G0  TO  819 
TYPE  5  —  CUMULUS  HUMULUS 
CTYPE(Kl)-5 

GO  TO  820 

TYPE  3  —  FOG/ STRATUS 

819  CTYPE(Kl)-3 
320  RETURN 


o  n  o  nno  n  n  n  n 


SUBROUTINE  PRECIP 

THIS  MODULE  IDENTIFIES  PRECIPITATION  INTENSITY  AREAS 

I- LIGHT,  2-MODERATE,  3-HEAVY,  AND  O-NO  RAIN. 

COMMON/ A/ ML( 128), CIRRUS, K.K1.ISEL 
COMMON/B/ CTYPE(64) ,PRCIPI(64) ,MSAVG<64) ,M 
REAL  ML , HE IGHT , CTEMP , MS AVG 
INTEGER  K, CIRRUS, PRCIPI ,CTYPE ,SUM1 

THESE  ARE  LIL.IAS  CONVERTED  GOES  THRESHOLDS 

II- 184 

12- 195 

13- 224 

ADD  IR  (ML)  AND  AVERAGE  VISUAL  (MSAVG) ,  THIS  SUM  WILL  BE  TESTED 
FOR  INTENSITY  OF  RAIN. 


SUM1-ML(K)+MSAVG(K1 ) 

IF(SUM1 .GE.I1 .AND.SUM1 ,LE.I2)GO  TO  43 
IF(SUM1 .GT. I2.AND.SUM1 .LE. I3)GO  TO  44 
IF(SUMl.GT.I3)GO  TO  45 
PRCIPI(Kl)-0 
GO  TO  49 

43  PRCIPIOCI)-! 

GO  TO  49 

44  PRCIPI(Kl)-2 
GO  TO  49 

45  PRCIPI(Kl)-3 
49  RETURN 

END 


SUBROUTINE  TEMHT 


C 

C  THIS  SUBROUTINE  CALCULATES  THE  CLOUD  TOP  TEMPERATURE  (CTEMP) 

C  USING  THE  IR  RADIANCE  FOR  THE  SURFACE  ( NCLR) ,  THE  AMOUNT  OF  CLOUD 

C  (ACDAVG),  EMISSIVITY  (E),  AND  THE  MEASURED  IR  RADIANCE  FROM  THE 

C  SATELLITE  (ML).  CTEMP  IS  THEN  CORRELATED  TO  AN  AREA 

C  REPRESENTATIVE  UPPER  AIR  SOUNDING  (UALV)  TO  DETERMINE  THE  CLOUD 
C  TOP  HEIGHT  (HEIGHT). 

C 

COMMON/ A/ML( 128), CIRRUS, K,K1,ISEL 

COMMON/ D/ HE IGHT ( 64 ) , CTEMP (  64  ) , ACDAVG ( 64 ) , E 1 , E  2 , ALAT , ALON , J i , 

*TLV( 11) ,E3 
COMMON/G/STOP(ll) 

REAL  ML , HEIGHT , CTE*fP , E 1 , E2 , MCLR.NCLD , ACDAVG ,  E ,  SIG 1 , SIG2 , SIG3 , 
*TLV,STDP 

INTEGER  I, CTYPE.PRCIPI, CIRRUS 

DATA  STDP/ 1010. 0,1000. 0,350.0,700. 0,500. 0,400. 0,300. 0,250. 0,200.0, 
*150.0,100.0/ 

E-0.0 

C  CHECK  TO  SEE  IF' THERE  ARE  ANY  CLOUDS,  IF  NOT  SKIP  TO  END 
IF(ACDAVG(K1 ) .LT.O. 1)G0  TO  61 
C  CONVERT  SURFACE  TEMPERATURE  (TSFC)  TO  NCLR 
IF(TLV( 1 ) ,GE.242.0)NCLR-2.0*( 331 .0-TLV( 1) ) 

I F( TLV( 1 ) . L  T . 242.0) NCLR-420 . 0-TLV( 1 ) 

C  IS  CIRRUS  FLAG  SET?  IF  YES,  USE  E2  FOR  E.  IF  NOT,  USE  El  FOR  E. 
I?( CIRRUS. SO.DGO  TO  850 

C  IS  THE  CD  OR  NS  CLOUD  FLAG  SET?  IF  YES  USE  E3  FOR  El. 

IF( 53.EO. 1 .0)00  TO  849 
E-El 

GO  TO  851 

849  E-E3 

GO  TO  851 

850  E-E2 

C  COMPUTE  CLOUD  RADIANCE 

851  NCLD-ML(K)-NCLR/(E*ACDAVG(K1))+>TCLR 

C  CONVERT  TO  CTEMP 

IF(NCLD.LE.178.0)CTEMP(K1)«331.0-NCLD*0.5 
I F( MOLD . GT . 1 78 . 0 ) CTEMP( K1 ) -420 . 0-NCLD 
’•TRITE  ( 12 ,353)NCLR,NCLD 
853  F0RMAT(/1X,5(1H*,E15.3),"NCLR,NCLD") 
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IF( ISEL.EQ.2)C0  TO  60 


INTERPOLATE  UPPER  AIR  SOUNDING  TO  CONTINUOUS  CURVE  OF  TEMP 
VERSUS  HEIGHT,  THEN  COMPARE  TO  CTEMP  TILL  GET  MATCH.  CHECK 
FROM  TOP  DOWN.  THEN  USE  MATCH  TO  SELECT  CORRECT  HEIGHT  FOR 
CLOUD  TOP. 

DO  812  18-2,11 
IA- 13-13 

IF(IA.EQ.ll. AND. TLV(11).GT.CTEMP(K1))G0  TO  316 
IF(TLV( IA) .LT. CTEMP (K1 ) )G0  TO  313 
IF(TLV(IA).EQ.CTEMP(K1))G0  TO  814 
312  CONTINUE 

813  IF(IA.EQ.2)IA-3 

RHG»EXP( ALOG( STDP( IA+1 ) ) -( ALOC( ST0p( IA+1 ) )-ALOG( STDP( IA) ) ) * 
*( TLV( IA+1 )-CTEMP(Kl ) )/( TLV( IA+1 )-TLV( IA) ) ) 

CO  TO  815 

816  TYPE  "TEMPERATURE  IS  ABOVE  100  MB" 

ACCEPT  "CONTINUE?", ICONT 
IF( ICONT.EQ.O) ISEL-8 

814  RHG-STDP(IA) 

815  HEIGHT(Kl)-RHC+0.5 

60  RETURN 

61  CTEMP (Xl)-TLV(l) 

HEIGHT(k.i)-STDP(1) 

RETURN 
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